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In this paper, we show how to discretize the abelian Chern-Simons gauge theory on generic 
planar lattices/graphs (with or without translational symmetries) embedded in arbitrary 2D closed 
orientable manifolds. We find that, as long as a one-to-one correspondence between vertices and 
faces can be defined on the graph such that each face is paired up with a neighboring vertex (and vice 
versa), a discretized abelian Chern-Simons theory can be constructed consistently. We further verify 
that all the essential properties of the Chern-Simons gauge theory are preserved in the discretized 
setup. In addition, we find that the existence of such a one-to-one correspondence is not only a 
sufficient condition for discretizing a Chern-Simons gauge theory but, for the discretized theory to be 
nonsingular and to preserve some key properties of the topological field theory, this correspondence 
is also a necessary one. A specific example will then be provided, in which we discretize the abelian 
Chern-Simons gauge theory on a tetrahedron. 


I. INTRODUCTION 

As the prototypical topological field theory, the Chern- 
Simons gauge theory, has had a deep and broad impact 
on a wide range of physics research, ranging from knot 
theorjffl and parity anomalies in quantum field theorjPl 
to the theory of the integer and fractional quantum Hall 
effect^® and the effective field theory description of chi¬ 
ral spin liquid^ in condensed matter physics (for a 
review see, e.g. Ref. 0 )- Although well understood 
as a continuum field theory, there is still limited under¬ 
standing on how to discretize this topological field the¬ 
ory on 2D lattices or graphs P This task turns out to 
be highly nontrivial. In particular, the topological and 
gauge-theoretic nature of the Chern-Simons gauge theory 
enforces strong constraints on the dynamics of the gauge 
fields. These constrains, if not treated carefully, can re¬ 
sult in inconsistencie s in t he discretized theory, making 
the theory ill-defined! 10 * 11 ! Until recently, the discretiza¬ 
tion has only been done only for a very special case, i.e., 
on a square lattice (with only nearest bonds) embedded 
in a torus. 10 11 It remains highly unclear whether simi¬ 
lar construction can be extended for other lattices, or for 
lattices embedded on other 2D manifolds aside from the 
torus, or in any discretized systems without translation 
symmetries (e.g. a graph). In a recent publication)- 12 ^ we 
presented a consistent construction of the Chern-Simons 
gauge theory on one of the simplest non-bipartite lat¬ 
tices in two-dimensions, the Kagome lattice, and used 
it to study the magnetizations plateaus of the spin-1/2 
frustrated quantum Heisenberg antiferromagnet on this 
lattice. 

The main purpose of the present paper is to develop 
a consistent discretization of the (abelian) Chern-Simons 
abelian U(l) gauge theory on general planar lattices and 
graphs. There are several motivations to search for a 
discretized Chern-Simons gauge theory on generic lat¬ 
tices/graphs. For example, it has been known that the 
Chern-Simons gauge theory plays a crucial role in the 


study of chiral spin liquid. Such an exotic state of matter 
can only be stabilized in the presence of strong geometric 
frustration. Much of the work in frustrated antiferromag¬ 
net s uses the fact that these systems are equivalent to 
a system of (generally interacting) hard-core bosons on 
the same lattice. The hard-core bosons are then mapped 
into a system of fermions coupled to a discretized Chern- 
Simons gauge field. ^Except for some very special exactly 
solvable models, in the study of such frustrated systems 
the dynamics and quantum fluctuations of the effective 
gauge fields are typically ignored, and frustrated quan¬ 
tum antiferromagnets are frequently described only at 
the level of the average field approximation! 10 14 How¬ 
ever, such a approximation is unreliable, and has a strong 
and obvious bias towards time-reveal breaking ground 
states. As shown in Ref. to correctly address the 
competition between different quantum ground states, it 
is necessary to go beyond the average field approximation 
by carefully introducing the correct quantum dynamics. 

It was recently realized that, in addition to the well 
known case of two-dimensional electron gases (in the con¬ 
tinuum) in strong magnetic fields, the fractional quan¬ 
tum Hall effect can also been stabilized in lattices even 
with zero net magnetic fieldP^HUl Fractional quantum 
Hall stat es ha ve also been explored on lattice systems 
earlier oi J 22 * 23 * and, more recently, also in optical lat¬ 
tice systemsp3K 27 lThis type of (discrete) fractional topo¬ 
logical states is now often referred to as the fractional 
Chern insulators or the fractional quantum anomalous 
Hall state. In particular, it has been shown that these 
(discrete) fractional Chern insulators are adiabatically 
connected to the corresponding fractional quantum Hall 
states in the continuum PsHSSl 

The Chern-Simons gauge theory is known to be the 
low-energy, hydrodynamic, theory of topological phases 
such as the fractional quantum Hall fluids,^ they are 
also expected to describe the low-energy and long¬ 
distance limit of topological Chern insulators, fractional 
or not EIHlIl Although Chern-Simons theories yield a nat- 
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ural description of the hydrodynamic behavior of topo¬ 
logical phases, apriori they are not absolutely neces¬ 
sary in the microscopic construction of a theory of such 
statesP However, all known fractionalized phases are in 
inherently strongly coupled systems and, notably in the 
case of the fractional quantum Hall fluids, the use of 
Chern-Simons gauge the ory in the microscopic deriva¬ 
tions has been i nvalu able.^ 4 Aside from some recent and 
promising workj^I 37 it is not yet clear what role does 
Chern-Simons gauge theory play in the theory of frac¬ 
tional topological Chern insulators. Although adiabatic 
continuity strongly implies that the theory of factional 
Chern insulators should be smoothly related to the the¬ 
ory of the fractional Hall effect on lattice systemsJ22123] 
where discrete Chern-Simons gauge theory is expected 
to be generally applicable. The general answer to these 
questions has remained problematic in view of the fact 
that the Chern-Simons gauge theory has not been dis¬ 
cretized on most of the lattices on which lattice fractional 
Chern insulators are known to occur (e.g. the checker¬ 
board lattice, the Kagorne lattice and the multiple-orbital 
square lattice). Although we will not give an answer for 
systems on lattices as general as it is needed, we will give 
an explicit construction for a large class of lattices, which 
include some of clear physical interest. 

In this paper, we propose and study a discretized 
Chern-Simons gauge theory on generic planar graphs em¬ 
bedded in arbitrary 2D closed and orientable manifolds.^ 
Same as in a lattice gauge theory, we will define the space 
components of the gauge field to live on the nearest- 
neighbor bonds of the lattice and the time components 
on the sites (vertices) of the lattice. We will consider 
only planar lattices (and hence with only non-crossing 
bonds).This is a lattice gauge theorj®^ but 0 ne with a 
broken time reversal and parity invariance. Here we will 
be interested in a version of Chern-Simons theory on a 
a class of spatial lattices with continuous time. Earlier 
work focused on the square lattice,®- 4414 ^ and recently we 
discussed the case of the Kagorne lattice? 42 Discretized 
versions of Chern-Simons gauge theory have been dis¬ 
cussed both in Euclidean space-time lattices,^ which suf¬ 
fer from the species doubling problem analogous to those 
of lattice fermions. By enlarging the scope of investi¬ 
gation from periodic lattices to graphs (with or without 
translational symmetries), our conclusions are generally 
applicable for a wide range of systems. 

We will require the discretized theory to retain the 
central features a topological field theory. Chern-Simons 
gauge theory on a continuous space-time manifold has 
several key features! 4 - It is gauge theory which means 
that it has a local symmetry under local (in space-time) 
gauge transformations. At the quantum level this re¬ 
quires that the quantum states of the physical Hilbert 
space be gauge invariant,^ and hence that the gener¬ 
ators of local time-independent gauge transformations 
must generate superselection sectors, i.e. the Gauss law 
is satisfied as constraint on the physical space of states. 
For this requirement to be consistently implemented, the 


generators of local gauge transformations must commute 
with each other on different spatial locations. In the case 
of the Chern-Simons theory this implies that the local 
magnetic flux must commute with each other (since they 
are the generator of time-independent gauge transforma¬ 
tions). This condition imposes stringent constraints on 
the possible form of the discretized theories,^ 114 and it 
is the main focus of this work. 

On the other hand, at the classical level, the Chern- 
Simons theory is topological in the sense that the action 
is invariant under general coordinate transformations and 
hence it is independent of the metric. A consequence of 
this feature is that the energy-momentum tensor classi¬ 
cally vanishes and, consequently, the Hamiltonian is also 
zero. Clearly, any lattice discretization implies a choice 
of coordinates, i.e. a fixed spatial metric. Furthermore, 
a change of the lattice stricture leads to a change in the 
form of the metric. Therefore, a lattice version of Chern- 
Simons theory cannot be explicitly independent of the 
metric and, in this sense, it cannot be formally topologi¬ 
cal. However, we will show below, that one can construct 
a Chern-Simons gauge theory for a large class of lattices 
a U(l) lattice gauge theory with continuous time (i.e. 
in “Hamiltonian” forn:P3) which is gauge-invariant. We 
will see that the resulting discretized theory neverthe¬ 
less has a vanishing Hamiltonian since the content of the 
action reduces to a set of (reasonably local) equal-time 
commutation relations and a set of local and commuting 
constraints. This theory is topological in the sense that it 
does not have local excitations, and that only the global 
degrees of freedom (non-trivial Wilson loops) matter. In 
the long-wavelength limit the discretized theory becomes 
(formally) the continuum Chern-Simons gauge theory. 

We find that such a Chern-Simons gauge theory can be 
constructed for arbitrary 2D planar graphs (lattices) as 
long as a local vertex-face correspondence can be defined 
on the graph/lattice. We adopt the following definition: 

Definition. A local vertex-face correspondence is a one- 
to-one correspondence between faces and vertices defined 
on a graph such that every vertex is adjacent to its cor¬ 
responding face (and vice versa). 

An example of such a correspondence is shown in Fig. |T| 

The relevance of this correspondence to Chern-Simons 
theory lies in the nature of this gauge theory. In the 
continuum the (abelian) Chern-Simons Lagrangian of a 
gauge field in 2+1 dimensions is (coupled to a matter 
current J^) 

Ccs[A] = ^e^ x A^Ax - (1.1) 

47T 

The Chern-Simons (CS) gauge theory is a topological 
field theory! 4 At the classical level the CS action is in¬ 
dependent on the metric of the manifold on which it is 
defined. The content of this Lagrangian, EqjOJ , is seen 
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FIG. 1. Part of a planar graph, on which a local vertex-face 
correspondence is defined. The disks and solid lines represent 
vertices and edges of the graph respectively. Each face is 
marked by a cross. The local vertex-face correspondence is 
indicated by dotted lines, that pair up each face with one (and 
only one) adjacent vertex. 


in Cartesian components 

£cs[-A\ — J 0 A q 

— —eijAidtAj — J ■ A (1-2) 

47T 

At the quantum level, the first term of the r.h.s. be¬ 
comes the requirement that the states in the physical 
Hilbert space, {|Phys)}, obey the “Gauss law” as a local 
constraint. Thus, the physical states are gauge-invariant 
and are annihilated by the generator of local gauge trans¬ 
formations, 


- k 

.27r 


B{x) - J 0 (x) 


jPhys) = 0 


(1.3) 


Hence, the physical states are those in which the local 
charge density Jq to the local magnetic flux B = eijdiAj 
are precisely related, i.e. flux attachment. The second 
term of the r.h.s. of Eq.( 1.2) implies that the components 
of the gauge field obey the equal-time commutation re¬ 
lations, 


[Ai{x),Aj{y)\ = i—eij6(x - y) (1.4) 

Further, the Hamiltonian of this system is zero unless 
sources are present, i.e. 

H = J-A (1.5) 

which is a consequence that a topological field theory 
does not have any excited states with finite energy. 

As with any lattice gauge theory, in a discretized 
Chern-Simons gauge theory the gauge fields (which are 
connections and hence are 1-forms) are naturally defined 
on the links of the lattice whil e the matter fields are de¬ 
fined on the sites of the latticeP®^ The field strength is 


a 2-form and it is defined on the elementary plaquettes of 
the lattice. While in a conventional lattice gauge theory 
the lattice is generally hypercubic (i.e. square in 2D), 
here we will consider more general (and translationally 
invariant) planar lattices. For instance, in Ref.j^j (and 
in Refs. .10] and HB) the Chern-Simons theory was de¬ 
fined on a square lattice and in Ref. [12] the theory was 
defined on a K agorne lattice. In both cases the Gauss 
law of Eq. (1.3) is naturally implemented as a constraint 
that relates the occupation number of a site (or vertex) 
to the gauge flux through a (uniquely defined) adjoint 
plaquette (or face). While in the case of the square lat¬ 
tice all plaquettes are identical (squares), in the case of 
the Kagorne lattice has three inequivalent sites in its unit 
cell and, correspondingly, three faces (two triangles and 
a hexagon) in its unit cell. Nevertheless, the correspon¬ 
dence of vertices to faces is one-to-one in both lattices. 


We will see here that this correspondence is a key fea¬ 
ture which will allow us to impose the constraint (and 
hence gauge invariance) in a unique way which, in ad¬ 
dition, does not break the point group (or space group) 
symmetries of the lattice. Below we will find a construc¬ 
tion of the Chern-Simons gauge theory on lattices for 
which for a charge located at a vertex, the magnetic field 
attached to it by the Chern-Simons gauge theory is lo¬ 
cated at the face that is naturally paired up with this 
vertex. 


Whether or not a local vertex-face correspondence can 
be defined for a graph is fully determined by the con¬ 
nectivity of the graph. In Sec. [TT[ we will provide an 
sufficient and necessary condition, which can be used to 
decide whether such a correspondence exists or not for 
an arbitrary graph. In Fig. [2} we show some examples of 
lattices that support such a correspondence (i.e. a dis¬ 
cretized Chern-Simons gauge theory can be constructed 
on these lattices). These examples include some of the 
lattices used in the study of chiral spin liquids and the 
lattice fractional quantum Hall effect (e.g. the Kagorne 
lattice). 

It is also worthwhile to emphasized that in the con¬ 
tinuum, the Chern-Simons gauge theory can be defined 
on arbitrary 2D closed and orientable manifolds. This 
plays a critical role in the phenomenon of topological de¬ 
generacy in fractional quantum Hall systems^. In addi¬ 
tion, it is also known that all the essential physics of the 
Chern-Simons gauge theory (in the continuum) is stable 
against the explicit breaking of the translational sym¬ 
metry, which is the underlying reason for the stability 
of the quantum Hall states against weak disorders. On 
the discretized side, however, it is still unclear whether 
the Chern-Simons gauge theory can be defined on any 
manifold aside from a torus and/or on a discrete graph 
without translational symmetries. Our study will provide 
an answer to these questions. 

In addition to those geometric considerations, a key 
consistency requirement of the gauge theory i s th at the 
lattice version of the local constraints of Eq. (1.3) must 
commute with each other and hence act as superselec- 
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tion rules on the Hilbert spacd^M (otherwise, these con¬ 
straints cannot be simultaneously satisfied). This consis¬ 
tency condition places restrictions on the commutation 
relations satisfied by the gauge fields defined on the links. 
For the square lattice this problem was solved by Eliezer 
and Semenoff,i° 11 and was more recently generalized by 
us to the case of the Kagome lattice.® In this paper we 
will show that the commutation relations can be defined 
consistently on any lattice (and graph) which obeys the 
one-to-one correspondence between vertices and faces. 
We will show that this restriction is implemented in terms 
of a suitably defined non-singular (and hence invertible) 
matrix. Therefore, the lattice Chern-Simons theory can 
be defined as a consistent gauge theory at the quantum 
level on these planar lattices and graphs. 

This paper is organized in the following way. In Sec. [TTJ 
we present a necessary and sufficient criterion for deter¬ 
mining whether a local vertex-face correspondence can 
be defined for an arbitrary graph/lattice. In Sec. Ill we 


write down the action of the discretized Chern-Simons 
gauge theory for generic graphs with a local vertex-face 
correspondence. In Secs. |IV||IX[ we prove that our dis¬ 
cretized gauge theory preserves all key features of the 
Chern-Simons gauge theory, including the gauge invari¬ 
ance, flux attachment, commutation relations, duality 
transformation and the locality condition. In Sec. [Xj 
we show that the existence of a local vertex-face corre¬ 
spondence is the necessary condition for discretizing the 
Chern-Simon theory, if we want the theory to be nonsin¬ 
gular and to preserve the key properties of the Chern- 
Simons gauge theory. In Sec. |XI| we present a simple ex¬ 
ample by discretizing the Chern-Simons gauge theory on 
a tetrahedron, which is a 2D planar graph on a sphere. 
In Sec. m we conclude our paper by discussing open 
problems and applications of this theory to a number of 
systems of interest. Details of the calculations are pre¬ 
sented in several appendices. 


II. THE LOCAL VERTEX-FACE 
CORRESPONDENCE 


We start our discussion by presenting all constrains 
and assumptions that will be enforced on the graphs 
(and lattices) that we will consider. In this paper, we 
study generic planar simple graphs embedded on arbi¬ 
trary closed and orientable 2D manifolds. Here, “planar” 
indicates that the graph can be drawn on a 2D manifold 
without any crossing bonds, while “simple” means no 
multiple bonds connecting the same pair of sites and no 
bond connecting a site to itself (See Fig. [TT] in Appendix 
[A] for an explicit example). The “simple” condition is au¬ 
tomatically implied for most (if not all) lattices studied 
in physics, while the “planar” condition holds for many 
(but not all) of them. 

For a planar graph G, we can construct the dual graph 
G* by mapping vertices to faces, and vice versa. Because, 
as will be discussed below, the dual graph will be needed 



(c) (d) 


FIG. 2. (Color online) Examples of lattices and graphs that 
support local vertex-face correspondences, (a) a square lattice 
with 1 vertex and 1 face per unit cell, (b) a Kagome lattice 
with 3 vertices and 3 faces per unit cell, (c) a dice lattice 
with 3 vertices and 3 faces per unit cell and (d) a lattice 
that contains 9 vertices 18 edges and 9 faces per unit cell. 
The (red) parallelogram marks a unit cell with lattice vectors 
indicated by the two (red) arrows. It is easy to verify that 
for all these lattices N v = Nf and for any subgraphs the 
number of faces never exceeds the number of vertices, which 
is a sufficient and necessary condition for the existence of (at 
least) one local vertex-face correspondence. 


for the dual gauge theory, we will also require the dual 
graph G* to be simple. For the original graph G, this 
condition implies that G cannot contain any dangling 
bonds, and that two faces in G can share at most one 
common edge. 

From now on, we will focus our study on graphs, on 
which a local vertex-face correspondence can be defined. 
Below, in Sec.|Xj we will prove that this constrain is nec¬ 
essary in order to preserve certain key defining properties 
of the Chern-Simons gauge theory. 

To determine whether a graph can support such a lo¬ 
cal vertex-face correspondence, we will use the following 
criterion: 

Criterion. A local vertex-face correspondence can be de¬ 
fined on a 2D planar graph G, if and only if the graph has 
the same number of vertices and faces (i. e. N v = Nf), 
and that for any subgraph of G the number of faces never 
exceeds the number of vertices (i.e. N' v > Nf). 

That this criterion is a sufficient and necessary condi¬ 
tion is proved in Appendix [B] The proof utilizes Hall’s 
marriage theorem by mapping the local vertex-face cor¬ 
respondence to Hall’s marriage problemW^ The marriage 
theorem is named after the British mathematician, Philip 
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(a )N V < N f 


(b)N v > N f 



(c )N f = N v 


FIG. 3. (Color online) Examples of lattices/graphs that do 
not support a local vertex-face correspondence, (a) a tri¬ 
angular lattice, which has 1 vertex and 2 faces per unit cell 
(N v < Nf), (b) a honeycomb lattice, which has 2 vertices and 
1 face per unit cell (N v > Nf) and (c) a lattice with N v = Nf 
but some of the subgraph has more faces than vertices, e.g. 
the dark area, which has 18 faces and 16 vertices. The (red) 
parallelogram marks a unit cell with lattice vectors indicated 
by the two (red) arrows. Each unit cell of this lattice contains 
27 vertices, 54 edges and 27 faces (N v = Nf). 


Hall, who should not be confused with the physicist Ed¬ 
win Hall, after whom the Hall effect is named. 

Using this criterion it is straightforward to determine 
whether or not a graph or lattice can support a local 
vertex-face correspondence. In Fig. 1 (Fig- [3]), we pro¬ 
vide examples of lattices/graphs, on which such a local 
correspondence exists (does not exist). In Fig. [3j the 
first two lattices do not support any one-to-one corre¬ 
spondence between vertices and faces, because the num¬ 
ber of faces does not match the number of vertices. The 
third example, shown in Fig. |3(c)| has the same number 
of faces and vertices and thus, in principle, a one-to-one 
correspondence between vertices and faces could be de¬ 
fined. However, in this case such a correspondence cannot 
be local, as proven in Appendix [B] because this lattice 
contains some subgraph, whose number of faces exceeds 
the number of vertices. For example, the dark area in 
Fig. |3(c)| shows a subgraph with 18 faces and 16 vertices. 

In Figs. |2(b)] and |2(c) the two lattices are dual to each 
other. Generically, if a graph G has a local vertex-face 
correspondence, so does its dual graph G*. This is be¬ 
cause one can construct such a correspondence for G* 
by simply swapping the vertices and faces in the origi¬ 
nal vertex-face correspondence defined on G. As a re¬ 
sult, our discretized Chern-Simons gauge theory always 
arises in pairs (one on the graph G and the other on the 


dual graph G*). In Sec. VI we will prove that these 
two gauge theories are dual to each other. This dual¬ 
ity relation is different from the continuum, in which the 
Chern-Simons theory is self-dual. A discretized Chern- 
Simons gauge theory is in general not self-dual, unless 
the underlying graph is self-dual. One example of a self 
dual graph is shown in Fig. |2(a)[ i.e., a square lattice. 
Another self-dual graph will be presented in Sec. |XI[ i.e., 
a tetrahedron. 

We conclude this section by highlighting some conven¬ 
tions adopted in this paper. For a graph G, we label the 
numbers of vertices, edges, and faces as N v , N e and Nf 
respectively, and we use the subindices v, e and, / to 
label each vertex, edge and face, respectively, where v, 
e and / take integer values (1 < v < N v , 1 < e < N e 
and 1 < / < Nf). For the dual graph G*, we will use 
the symbol to label every object. For example, ver¬ 
tices, edges, and faces of the dual graph are labeled as 
v*, e* and /*, respectively. In addition, for convenience, 
if a vertex v in graph G is mapped to the face f* in the 
dual graph, we will use the same integer to label them, 
i.e., v = /*. Same is true for corresponding e and e* (/ 
and v*). Throughout the paper, repeated indices will be 
summed over unless specified otherwise. For the gauge 
field, the time-component lives on vertices and thus will 
be labeled as A v . The spatial components (i.e. the vector 
potential) are defined on edges, and thus will be shown 
as A e . Because the vector potential is a vector, we must 
choose a positive direction for each edge (from one of its 
end to the other). The vector potential A e on an edge e 
is positive (negative), if it is along (against) the direction 
of the edge e. In graph theory, after a direction is as¬ 
signed to each edge, the graph is called a directed graph 
(or a digraph) S3 


III. THE DISCRETIZED CHERN-SIMONS 
ACTION 

In this section, we construct the action of the dis¬ 
cretized Chern-Simons gauge theory. We should empha¬ 
size that as long as the conditions discussed in the pre¬ 
vious section are satisfied, our construction is applicable 
for arbitrary graphs. In addition, as will be shown below, 
the discretized gauge theory obtained here is a topolog¬ 
ical field theory, whose action only relies on the connec¬ 
tivity of the graph without any free parameter, except 
for a quantized topological index k. 


A. The M matrix and the K matrix 

In this section, we define two matrices for arbitrary 
graphs with a local vertex-face correspondence. For a 
graph satisfying the criterion given in the previous sec¬ 
tion typically there is more than one way to define the lo¬ 
cal vertex-face correspondence, and different choices here 
will in general result in different M and K matrices and 
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thus lead to slightly different actions. Here we choose an 
specific (albeit arbitrary) one, consistently throughout 
the lattice. 

The vertex-face correspondence defines a matrix M v j 
with dimensions N v x Nf. The first index of this matrix 
runs over all vertices, while the second one indicates faces 
in the graph. If vertex v and face / are paired-up accord¬ 
ing to the vertex-face correspondence, then M v j = 1. 
Otherwise, the matrix element is zero. Hence, 

/ 1 if v is paired with / , > 

(0 otherwise ' ' ' 

Because the vertex-face correspondence requires N v = 
Nf. the matrix M is a square matrix. In addition, it is 
easy to realize that, by definition, M is an invertible and 
orthogonal matrix, i.e. the inverse matrix M _1 is the 
transpose matrix, M T = Af _1 . 




(c)??i = — 1 and 772 = +1 (d)? 7 i = — 1 and r ]2 = —1 


which we will denote by K , whose two indices run over all 
edges of the graph (with N e being the number of edges), 

j. _ J ±| if e and e' belongs to the same face 
e,e ' | 0 otherwise 

(3.2) 


If there exists a face / such that e and e! are both edges 
of this face, the component of the matrix K e e f is ±1/2. 
Otherwise the matrix element vanishes. For nonzero 
K ee >, the ± sign is determined by the following formula, 


v r h x r) 2 1 

1\ ( , P ! =-= zb 


2 ’ 


(3.3) 


where 77 1 = ±1 and 772 = ±1 are two Z 2 integers. 

The sign of 771 is determined using the following rule. 
As shown in Fig. [4j we first mark the vertex that is paired 
up with / in the local vertex-face correspondence using 
a (red) circle. After that, we go from the edge e to the 
edge e 1 by moving counter-clockwise around the face /. 
If the path goes though the specially marked vertex (the 
red circle in Fig. [4]), 771 = ±1, and otherwise 771 = —1 

The sign of 772 is determined by the directions of the 
two edges e and e!. As discussed above, to define the 
vector potential, we must specify the direction for each 
edge. When we goes around the face / in the counter¬ 
clockwise direction, if both e and e! are pointing along (or 
opposite to) the direction of the path, 772 = ±1. If one of 
them points along the path while the other is opposite, 
V 2 = - 1 - 

With 771 and 772 , their product (multiplied b y -1) , —771 x 
772 = ±1, determines the sign of K e y in Eq. ( |3.2| ). Some 
examples can be found in Fig. [4j 


B. The Action 

With the two matrices defined above, we can now write 
down the action of our discretized Cliern-Simons gauge 
theory: 


FIG. 4. (Color online) Nonzero components of the K matrix. 
Here, we consider two edges e and e.', which belongs to the 
same face / (otherwise K ee i = 0). Based on the local vertex- 
face correspondence, the face / is paired up with one of its 
vertices, which is marked by the (red) circle. We go around 
the face / from e to e ' by following the direction of the positive 
orientation marked by the (blue) circle at the center of the 
face. In Fig. (a) and (b), the path from e to e' goes through 
the special site (marked by the red circle), and thus 771 = ± 1 . 
For Figs, (c) and (d), the special site is not on our path, and 
thus 771 = —1. The sign of 772 is determined by the orientation 
of e and e . If their directions are both along (or opposite) 
to the direction of the positive orientation [Figs (a) and (c)], 
772 = ±1. Otherwise [Figs, (b) and (d)], 772 = —1. Once 771 
and 772 are determined, the value of K ee i can be obtained as 
K e y = —771 x 772/2 = ± 1 / 2 . 

In addition to M, the local vertex-face correspondence 
can be used to define another N e x N e square matrix, 


S= N~ 

27r 


dt 


A v M v j<f>f — -A e K ee iA e 


(3.4) 


Here, we sum over all repeated indices. The index v, f 
and e run over all vertices, faces and edges respectively. 
A v is the time-component of the gauge field, which lives 
on vertices and A e represents the spatial components, 
which are defined on edges. Here, A represents the time 
derivative, K and M are the two matrices defined in the 
previous subsection, and is the magnetic flux on the 
face /, which equals to the loop integral of A e around /, 


4*/ — £f, e A e . (3-5) 


Here we sum over all edges and 


£/,e 


±1 e is an edge of / with positive orientation 
— 1 e is an edge of / with negative orientation 
0 e is not an edge of / 

(3.6) 
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The sign of £/, e is determined by going around the face 
/ along the counter-clockwise direction. If the direction 
of the edge e is along this path, £f e = +1. Otherwise, 
= —1. As can be seen from Eq. (3.5), the matrix 
is a discretized curl operator (Vx) for planar graphs. 

On a square lattice, the action that we constructed here 
reduces to the action found in Refs. [H and HU which 
can be considered as a special situation of our generic 
construction. Similarly, for the Kagome lattice this gen¬ 
eral construction reduces to the construction that we pre¬ 
sented in Ref. m- 

We conclude this section by comparing our discretized 
theory with the Chern-Simons gauge theory in the con¬ 
tinuum. For comparison, we choose to write down the 
action in the continuum in a special form 


S = 


2n 


dtdx 


AqB 2 AjCij'Aj 


(3.7) 


Here A 0 is the time component of the gauge field. A, 
and Aj are the spatial component with i and j being x 
or y. ej j is the Levi-Civita symbol and B is the magnetic 
held perpendicular to the 2D plane. The first term here 
enforces the flux attachment and the second term dictates 
the dynamics of the vector potential A x and A y . 

By comparing Eq. (3.4) with Eq. (3.7), we find that 


our discretized theory is in close analogy to the contin¬ 
uum case. Here, the M matrix dictates the flux attach¬ 
ment (i.e. Gauss’ law) and the AT-matrix plays the role 
of the Levi-Civita symbol. It is worthwhile to highlight 
that, just as the Levi-Civita symbol, the K matrix is 
antisymmetric 


K ,= —K / 

,e5 


(3.8) 


This can be verified easily by noticing that 771 — > —r)i 
and 772 —> 772 , if we swap e and e!. This antisymmetry 
property is in fact expected. If we look at the second term 
in our action, Eq (3.4), because J dtA e A e i = — f dtA e A e i 
(integration by part), only the antisymmetric part of K 
contributes to the action. 

In the next six sections, we will demonstrate that our 
action indeed offers a discretized Chern-Simons gauge 
theory on generic graphs by showing that all the key 
properties of the Chern-Simons gauge theory are pre¬ 
served by our action. 


IV. GAUGE INVARIANCE 

For a gauge theory, the action must be gauge in¬ 
variant. In this section, we will verify that our action 
[Eq. (3.4)] preserves the gauge symmetry. In the case 
of Chern-Simons, this is also true provided the mani¬ 
fold has no boundaries. Furthermore, invariance under 


large gauge transformations (which wind around non- 
contractible loops of the systems) holds only if the index 
k is an integer. 1 These conditions are satisfied by our 
discretized Chern-Simons theory. 


A. Gauge transformation on a graph 

For a graph/lattice, a gauge transformation takes the 
following form 


A v — >A V — dt<j 
A e — >A e — D v 


(4.1) 

(4.2) 


where <f> v is an arbitrary scalar function defined on ver¬ 
tices. The first formula [Eq. (4.1)] is the gauge transfor¬ 
mation for the time component of the gauge field, while 
the second one [Eq. (4.2)] is for the spatial components. 
The matrix D ve in Eq. (4.2) is the incident matrix of the 

grapfpl 


. +1 if v is the positive end of e 
Dy v = < —1 if v is the negative end of e 
0 otherwise 


(4.3) 


Here, we called the vertex v a positive (negative) end 
of the edge e, if v is one of the two ends of e and the 
direction of the edge e is pointing towards (away from) v. 
The incident matrix contains all the information about 
the connectivity of the graph, as well as the direction 
assigned to each edgeP® The incident matrix plays the 
role of a (discretized) gradient operator, V, which can 
be seeing easily by noticing that 


Dy 


(4.4) 


where </>„ is an arbitrary scalar function and the edge e 
points from v-i to v\. As a result, Eq. (4.2) can be con¬ 


sidered as a discretized version of A —> A — V0. Later, 
we will show that the incident matrix also serves as a 
discretized divergence, V-. 


B. Gauge symmetry 

As proven in Appendix [C] the sufficient and necessary 
condition for the action of Eq. ( |3.4[ ) to be gauge-invariant 
is that the following identity is satisfied 


A^u,/£/,e — D v ^ e 


(4.5) 


where £/ je is defined in Eq. (3.6) and the incident matrix 


D v e / is defined in Eq. (4.3). In this section, we prove that 


this condition is indeed valid for the M and K matrices 
constructed in Sec. MM 


To verify Eq. (4.51, we need to prove that the rela¬ 


tion holds for any e and v. Here, we classify all possible 
situations into three cases: 

1. e and v don’t belong to the same face. 

2. e and v belong to a same face but v is not an end 
of e 

3. v is an end of e 


Here we verify Eq. (4.51 for each of these three cases. 






























FIG. 5. Gauge invariance of our theory for (a) Case II and 
(b) Case III. In Fig. (a), we marked two additional edges of v, 
ei and e 2 , which are edges of /i. In Fig. (b), we labeled two 
faces /i and /2 and two additional edges ei and e 2 , such that 
e is the common edge shared by /i and / 2 , while ei and e 2 are 
two edges of v, which are adjacent to fi and / 2 respectively. 
Dashed lines represent (possible) additional edges of v, which 
are irrelevant for our proof and thus are not labeled. Although 
we assume a specific set of orientations for edges in these two 
figures, none of our final conclusions relies on the choice of 
orientations for each edge, as proven in Appendix [C| 


1. Case I 


The first case, where e and v don’t belong to the same 
face, represents the situation where e and v are separated 
far away from each other. It is easy to verify that in this 
case both sides of Eq. (4.51 vanish. 


For the l.h.s., M v j ^ 0 requires v being a vertex of the 
face / and £/ je ^ 0 requires e being an edge of /. For Case 
I, these two conditions cannot be satisfied simultaneously, 
and thus M v j£f, e = 0. 

For the r.h.s., K e e i ^ 0 implies that e and e' are edges 
of the same face, which will be called the face /. If 
D V}e ' 0, v must be one end of e ', which means that v 
is a vertex of /. As a result, to get a nonzero K e yD v y, 
both e and v must both belong to the same face /. This 
is in contradiction with the assumption of Case I, and 


Because both sides of the equation are zero, then 


Eq.( 4.5) holds for Case I. 


are labeled as e± and e 2 in Fig. |5(a)[ 

Using the edge directions shown in Fig. |5(a)| it is easy 
to verify that 


K — — 

1Y e,ei — 


K — — 
IY e,e2 


?7l;e,ei ??2;e,ei 


^7l;e,e2^2;e,ei 


m. 


2 

^l;e,e2 


Dy,ei - 1 


J-^v,e2 — 


1 


e ’ ei (4.6) 

(4.7) 

(4.8) 

(4.9) 


and thus 

K e ,e' D v ,e r — -^-e,ei H - -^e,e2^v,e2 — 


(4.10) 


Here, we shall distinguish two different situations: 1) 
v is paired up with /i according to the vertex-face corre¬ 
spondence, and 2) v is not paired up with /i. 

If v is paried up with /i, M v j vanishes for all /, except 
for f = fi, and therefore, 

Mi ;,/£/,e = M v ,f!^f u e = +1- (4-11) 


Here, we don’t sum over the repeated index /i and we 
used the fact that M v j 1 = 1. For the orientation shown 
in Fig. 5(a)| £/, e = +1. For the r.h.s. of Eq. (|4.5|), 
we can use Eq.(4.10). Here, we have ?7i ;e ,ei = +1 and 
r ?i;e,e 2 = — 1, and thus, 


K e , e ,D e ,, v = ?7l;e ’ ei 2 ??1;e ' e2 = \ + \ = +l. (4.12) 


By comparing the two equations above, we find that 

Mi >,/£/,e = fC,' 

If v is not pair ed up with / 1; M v j 1 £ : f i e = 0. For the 
r.h.s. ofEq. (4.5), it is easy to verify that =7 ?i;e,e 2 ) 

and thus 




Again, we verified Eq. (4.5). 


(4.13) 


2. Case II 


3. Case III 


The second case, where e and v belong to one same 
face but v is not an end of e, is shown in Fig. |5(a 3 In 
this figure, without loss of generality, we choose a spe¬ 
cific dire ction for each edge. As proved in Appendix |D| 
Eq. (4.5) is independent of the choice of the edge direc¬ 
tions. Therefore, although we only consider one specific 
direction arrangement here, the conclusion is generic. 


In Fig. |5(a)[ both v and e belong to the same face, f±. 
Because v is a vertex of the face /i, two of the edges of 
the face /i must have v as their end. These two edges 


For the last case, shown in Fig. |5(b)| because each 
edge in our graph is shared by two and only two faces 
(as shown above in Sec. 0 , we can label the two faces of 
the edge e as fi and / 2 . In addition, we also label two 
edges of v, e\ and e 2 , where e\ is an edge of fi and e 2 is 
an edge of / 2 . Because we have assumed that the graph 
and the dual graph are both simple (see Sec. 0, ei ^ e 2 . 

Same as in Case II, here too we only need to consider 
one specific set of directions for the edges and the con¬ 
clusion will be generic. Using the directions shown in 
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Fig. |5(b)| we have 



^7l;e,ei ^2;e,ei ^7l;e,ei 

(4.14) 

2 2 

K — — 
IY e,e2 — 

?7l;e,e2^72;e,ei ^7l;e,ei 

(4.15) 

2 2 

1 

1 

(4.16) 

Dy,e 2 — 

1 

(4.17) 

and therefore 



K / D 

1Y e,e J - y v y 

e' — -^e,ei-^y,ei -^e,e 2 -^v,e 2 



^7l;e,ei “1“ ^7l;e,e2 

2 

(4.18) 


Again, we distinguish two possible situations: 1) v is 
paired up with /i or fi according to the vertex-face cor¬ 
respondence 2) v is not paired up with either /i or / 2 . 

For the first situation, without loss of generality we 
assume that v is pa ired up with /i. Using the directions 
shown in Fig. |5(b)| we find that 

Mv,f£f,e = = Al- (4-19) 


Here, we don’t sum over the repeated index f\. In addi¬ 
tion, we also have ?7i ;e ,ei = 7 7i;e,e 2 = +1, and therefore 


Ke,e'D e ,, v = T?1;e ' ei + 7?1;e ’ 62 = +1. (4.20) 

So, we find = K e yD e > tV . 

If v is not paired up with either /i or / 2 , M v j^f^ e = 0, 
because it is impossible to make both M v j and 
nonzero. It is also easy to verify that here ?7i ;e ,ei = 
— ’ r ll;e,ei and thus 


K e e >D e r v = ??1;e ’ ei + ??1;e ' e2 = 0. 


(4.21) 


Once again, we get M„ j£/ >e = K e ^D e ^ v . 

By summarizing all p ossib le situations discussed 
above, we have verified Eq. (4.5). Therefore, we conclude 


that our theory is invariant under local gauge transfor¬ 
mations. 


V. FLUX ATTACHMENT 


states be invariant under local time-independent gauge 
transformations,® as we discussed in the Introduction, 
c.f. Eq.(1.3). This condition requires regularization (in 
the form of splitting the position of the charge and the 
flux) which leads to a pr oper fr aming of the knots rep¬ 
resented by Wilson loops. 1 ®-® 0 For a discrete system, 
however, because electric charges live on vertices, while 
magnetic fluxes are defined on faces (which takes care of 
the regularization), it is necessary to specify one addi¬ 
tional rule to dictate the location of the magnetic flux 
for charged particles at each site. This is achieved by 
the local vertex-face correspondence introduced in Sec. |T] 
Here too, this constraint amounts to the conditions that 
the states of the gauge theory be invariant under time- 
independent gauge transformations.® 

Because our action, Eq. (3.4), does not contain any 
dynamics for the time component of the gauge field A v 
(just as in any gauge theory), A v is not a dynamical 
field but its role is to enforce a constraint.® By taking 
a variational derivative of A v , we get the charge at the 
vertex v, 


qv ~5A v ~ 27t M "’ /$/ ’ (5 ' 2) 

which is proportional to the magnetic flux in the face /. 
Because M is an orthogonal matrix, this equation implies 
that 


2tt 

$/ — —q v M v j. 


(5.3) 


This equation is the discrete version of the flux attach¬ 


ment, analogous to Eq. (5.1). 


Here, we find that for a charge at a vertex v , a mag¬ 
netic flux is bound to it and the flux is located at the 
face /, which is the partner of v according to the vertex- 
face correspondence. This is the physical content of the 
vertex-face correspondence. 

We conclude this section by emphasizing that the flux 
attachment rule here is local , because we have required 
the vertex-face correspondence to be local , i.e., the mag¬ 
netic flux attached to a charge is located on a neighboring 
face. For a discrete system, this set up off ers the closest 
analogy to the delta function in Eq. (5.1). 


A key property of the Chern-Simons gauge theory is 
the constraint of flux attachment, which binds a magnetic 
flux with each charged particle. For a point charge q at 
location ro, the corresponding magnetic field is 

On t 

B{r) = ~k 9 ^ ^ - r °)’ ( 5 - 1 ) 


In the continuum classical theory, the flux and the charge 
are located at th e sa me position, as indicated by the 5- 
function in Eq. (5.1). In a continuum quantum gauge 
theory this condition is a constraint on the physical 
Hilbert space, and is the requirement that the quantum 


VI. DUAL GRAPH, DUAL THEORY AND THE 
INVERTIBILITY OF THE K MATRIX 


In this section, we verify two key (and essential) prop¬ 
erties of the discretized Chern-Simons gauge theory: 

1. The AT-matrix is invertible 

2. For any discretized Chern-Simons gauge theory 
constructed above, one can construct another dis¬ 
cretized Chern-Simons gauge theory on the dual 
graph. 
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Later, we will prove in Appendix [E] that the theory de¬ 
fined on the dual graph is in fact the dual theory of the 
original discretized Chern-Simons gauge theory. 

As has been addressed in literature, the AT matrix 
must be nonsingular (invertible) in order to ensure the 
correc t dyn amics for a discretized Chern-Simons gauge 
theory.^® 1 One way to realize this is by noticing that 
the inverse of the K matrix offers the commutator of the 
vector potential A e (see Sec. VII for more details), and 
therefore, to avoid singularities in the commutator, the 
K matrix must be invertible. 

Here, we will first verify the second property listed 
above by directly constructing a Chern-Simons gauge 
theory on the dual graph in Secs. |VI A| and |VIB[ Then, 
in Sec. |VI C[ we prove that K is invertible by finding di¬ 
rectly the inverse matrix of AT, which is in fact the K* 
matrix defined on the dual graph with a minus sign. Fi¬ 
nally, in Sec. |VI D| as a consistency check, we prove that 
the gauge invariance condition for the original graph and 
that of the dual graph are actually equivalent to each 
other. 


A. Duality transformation 


For a planar graph G, one can construct the dual graph 
G* by putting a vertex v* in each face of G and then 
connecting two vertices in G* if their corresponding faces 
in G share a common edge. It is easy to check that the 
dual of a dual graph is the original graph (G*)* = G. For 
the lattices shown in Fig. [2j the square lattice is self-dual, 
while the Kagome lattice and the dice lattice are dual to 
each other. 

For simplicity, we will use the same integer to label / 
and v*, if / is mapped to v* under the duality transfor¬ 
mation. Similarly, we use the same integer to label e and 
e* (v and /* ), if they are dual to each other. In addition, 
we choose the direction for each edge in the dual graph 
such that n e x n* > 0, where n e and n* are unit vectors 
along the direction of the edge e and its dual edge e*. In 
other words, we rotate the edge e counter-clockwise until 
it aligns with e*, and then the direction of the rotated 
edge e determines the direction of e*. 

With this convention, the incident matrix of the dual 
graph D *„ e . coincides with the £/ ie matrix of the original 
graph, Eq. (pL6|), 


D* 


= £/,e- 


( 6 . 1 ) 


Similarly, the e » matrix for the dual graph is in fact 
the incident matrix of the original graph D 1 up to an over 
all minus sign 


£/*,e* = ~Dv,e. ( 6 - 2 ) 

Here, we require v* = / and e* = e as shown in the 
previous graph. The physics meaning of these two rela¬ 
tions is that if the duality transformation maps a face / 
of a graph G into the vertex v* in the dual graph G*, 


then a loop around the face / is mapped to all the edges 
connected to vertex v*. and vice versa. 

It is easy to realize that under a duality transforma¬ 
tion, the local vertex-face correspondence in the origi¬ 
nal graph is transformed into a local vertex-face corre¬ 
spondence in the dual graph. As a result, we can use 
exactly the same construction to obtain a discretized 
Chern-Simons gauge theory on the dual graph 



dt 


- -A* e ,K* e , e „A* e 


(6.3) 


Here, A*„ and A*» are gauge fields defined on the dual 
graph with = £5, .A* being the magnetic flux of 
this gauge field on face /*. The M* and K* matrices 
are constr ucted using the same rules discussed above in 
Sec. Ill A In Appendix. |Ej we show that if k* = — 1/k , 
this action is the dual theory of the original discretized 
Chern-Simons gauge theory, Eq. ( |3.4[ ). 

It is straightforward to verify that the M* matrix is the 
transpose of the M matrix. Because M is an orthogonal 
matrix, it implies that M* is the inverse of M 


M* =M T = M~ 


(6.4) 


Below, we will study the K* matrix and prove that it is 
the inverse of the AT matrix up to an overall minus sign. 


B. the AT* matrix 


In this section, we show that the K* matrix can be con¬ 
structed directly in the original graph G, without going 
to the dual graph G*. This construction is fully equiva¬ 
lent to the dual-graph construction used in the previous 
section. However, as will be shown in the next section, 
by constructing the K* matrix in the original graph, it 
is more convenient to study the relation between the AT 
matrix and the K* matrix. 

As mentioned above, we label each edge in the dual 
graph using the same index as the corresponding edge 
in the original graph, (i.e. e* = e). Therefore, we can 
rewrite the AT* matrix using the edge indices of the orig¬ 
inal graph (e and e'), 


= K* e , (6.5) 

where e and e' are edges of original lattice and they are 
dual to e* and e'* respectively. The matrix AT* e , is now 
defined on the original graph, and thus we can translate 
the definition of the AT* matrix to the original graph. 
Using the original graph, it is straightforward to verify 
that 


AT* 


± j if e and e! share a vertex 
0 otherwise 


( 6 . 6 ) 


If e and e! do not share a common endpoint, K e ^ = 0. 
Otherwise, 


A e,e' = - 2- = ± 2' 


(6.7) 
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where 77 * = ±1 and 772 = ±1 are two Z 2 integers. To 
determine the sign of 77 *, we first label the common end 
of e and e! as v. Under the vertex-face correspondence 
(of the original graph), v is paired up with a neighboring 
face /. Now, we go from the edge e to the edge e' by 
moving around v in the counter-clockwise direction. If 
the path goes through the face /, 77 ^ = + 1 , and otherwise 
77 * = —1. The sign of 77 ^ is determined by the directions 
of edges e and e'. If both of them point toward (or away 
from) v , 77 ! = + 1 , and otherwise 772 = — 1 . 


C. I<* = -K- 1 


In this section, we prove that 

K* = —K~ l 


( 6 . 8 ) 


and thus K is invertible. To prove Eq. (6.8), we shall 
verify the following relations 


K e , e "K* e> , e , = K* 


,„K e 


= -5e 


(6.9) 


where 5 e ,e' is the Kronecker delta. In this section, we will 
only demonstrate K e e "K*„ e , = — S e y, while one can use 
the same method to prove K* e „K e "y = —S e y. 

Here, we need to consider six different cases, 



(a) Case 1 (b)Case 2 



1 . e = e'. 

2 . e / e', and e and e' share an endpoint, and e and 
e' are edges of the same face. 

3. e / e 7 , and e and e' share an endpoint, but e and 
e' are not edges of the same face. 

4. e and e' do not share any endpoint, but the belongs 
to the same face. 

5. e and e' do not belong to the same face, but there 
is a face /, where e is an edge of / and one of the 
endpoints of e' is a vertex of /. 

6 . otherwise 


Among all the six cases, 5 e ^ = 1 for the first one, and 
^e,e' = 0 for all others. In Fig. [gJ we show the first five 
cases. Here, we mark e, e' and all other edges that con¬ 
tribute to K e e n K*„ e , using solid lines. Other (possible) 
edges, which do not contribute to K ete nK*„ e ,, are la¬ 
beled as dashed lines. Using Fig. [ 6 j it is easy to notice 
that for all the first five situations 


K ,, K*„ 


= Y, K ^ K l,e' 

i 

= ^ ^ r h;e,eiV2;e,e i ' l ll;ei,e' r l2;ei,e l 
i 


( 6 . 10 ) 


FIG. 6. Possible cases for computing K ete iiK* r , e /. Solid lines 
marks the edges e, e! and all edges that will contribute to 
K, : e "K*n e ,. Dashed lines are (possible) additional edges, 
which don’t contribute to K f , f ,nK'^n e /. 


1. Case 1 


For the first situation [e = e! as shown in Fig. 6(a)] , 
we have 


K e ,e"Kl„ e , =e = Y, K e ,eJ C„e- (6-11) 


Here we don’t sum over repeated indices on the r.h.s. 
of the equation, a nd th e four edges ei, e 2 , e 3 and e 4 
are marked in Fig. | 6 (a)[ For any e^, we can verify that 
771 = 77 * and 772 = — 772 - Therefore, Ke,e t = — K*. e . As a 
result, 


JU r"K*„ e ,_ / , = Y (^e.ej 2 = - £ ( ±9 ) = -1 

( 6 . 12 ) 


i=1 


i =1 


and thus, we find that K e yK*„ e , = —5 e ,e' f° r e = e!. 


2. Case 2 


Here, for each e,, 771 , 772 , 77 ^ and 772 are obtained using _ 

the rules defined above. Below, we compute K e yiK*„ e , For the second case, shown in Fig. | 6 (b)[ it is straight- 
for each situation using this formula. forward to verify that if r] 2 - e ,ei an d 7 72 ;e,e 2 have the same 
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sign, then e\ and e-i must both point towa rds (or away 
from) the two vertices shown in Fig. |6(b)| As a result, 
^ 2 -ei e' an d 7 ?2-e 2 e' mus t have opposite sign, and thus 


r l2-e,eiV2-ei,e l r l2',e,e 2 V2;e 2 ,e' 


(6.13) 


Similarly, we can show that if ? 72 ;e,ei = ~V 2 ;e,e-? > we must 
have 772- ei e , = r? 2 ;e 2 e'- Ami therefore, Eq. (6.13) is aways 
valid for Case 2. 

For r)i s, we need to examine three different cases. Here, 
we consider the face / formed by e, e', e± and (possibly) 
other edges, and ask whether the vertex-face correspon¬ 
dence pairs up / with one of these two vertices. In gen¬ 
eral, there are three possibilities 


1. / is paired up with the vertex on the left; 


2. / is paired up with the vertex on the right; 


3. / is not paired up with either of them. 

For the first sitaution, we have ?7i ;e ,ei = — 1 , V*- ei e' = — 1 
and ?7i;e,e 2 = r li-e 2 e > ■ Therefore, we find 


r ll-,e,ei'ni-,ei,e' ~ r h',e,e 2 Vl;e 2 ,e' ■ 


(6.14) 


For the second situation, we have r/i ;e , ei = — 1, Vi- ei e’ = 


+1 and r]i- ete2 = — Vi-e 2 e' ■ Therefore, Eq. (|6.14[) is still 


valid. For the third situation, it can be shown that 
Vl;e, e± = +1, Vl- ei ,e' = + 1 and fo;e,e 2 = ? ?l;e 2 ,e'- Thus, 
Eq. (6.14| is still valid. 


In summary, we find that Eqs. (6.13) and (|6.14|) always 
hold for this case, Fig. |6(b)| By mu i 
equations together, we get 


tiplying these two 


4- Case 6 

The last case, 6, is easy to verify, because here e and 
e' are far away from each other, so that for any e", either 
K ei& n or K*„ e , is zero. Therefore, K e e n K*,, e , = 0 = 

$e,e' ■ 

By summarizing all the possible cases, we conclude 
that KK* = —I. We can use the same method to prove 
that K*K = — I and thus K* = —A' -1 . This result also 
proves that the K and K* matrices that we constructed 
above are invertible. 


D. Gauge invariance in the dual graph 


As shown above in Sec. |IV[ in the original graph, our 
action of Eq. (3.41 is gauge invariant if and only if 




( 6 . 20 ) 


For the dual graph, there is a similar condition for the 
gauge invariance. 


= K,^D* v%e „ 


( 6 . 21 ) 


In this section, we prove that these two conditions are in 
fact equivalent as long as M* = M -1 and K* = —1\ -1 


We start from Eq. (6.21) and change the dual graph 


(face, edge or vertex) labels into the corresponding labels 
of the original graph 


M f, v D v ,e = Af,e' 


( 6 . 22 ) 


K K* , = -K K* , 

ly e,e 1 * y e 1 ,e' e > e 2 ^ 1 e 2 ,e' 


(6.15) 


Utilizing Eq. (6.10), this relation implies that 
, = ( 

,=-5 e 


K ee nK*,, e , = 0, in agreement with the relation 

Ke,e"K e --,e 


and here we also us e th e rela tion s D*. e , = £/ i6 and 
= -D v>e , Eqs. |G~T| ) and ( |fo2| ). 

If M* = M _1 and K* = —AT” 1 , the formula above 
implies that 


— K e ,e'£f,e 


(6.23) 


3. Cases 3, 4 an d 5 


Using the same approach, we can show that for the 
third and the fourth cases, shown in Fig. 6(c) and |6(d)[ 


=%;e,e 2 ?72;e 2 ,e'» ( 6 - 16 ) 

T h;e,e 1 Vl- ei ,e' = ~ Vl;e,e 2 Vt,e 2 ,e' ■ ( 6 - 17 ) 

Once again, we get K e yK*„ e , = 0 = —5 e y. 

For the fifth case, shown in Fig. |6(e)| we have 


? ?2;e,ei%;ei,e' = - 5 ?2;e,e 2 ??2;e 2 ,e') ( 6 - 1 8) 

r h-,e,e 1 vt,e 1 ,e' =Vl-,e,e 2 V*l-e 2 ,e' ■ ( 6 -!9) 


Thus, K e ,e"K* e „ e , = 0 = — S e 


By multiplying the matrices M and K on both sides, we 
recover the condition of gauge invariance in the original 
graph, Eq. (4.5). Therefor e, we find t hat th e two gauge 
invariance conditions, Eq. (6.20) and (6.21), are equiva¬ 
lent. 


VII. COMMUTATION RELATIONS AND 
THE A " 1 MATRIX 


The Chern-Simons theory in the continuum has a very 
special commutation relations. In particular, the commu¬ 
tator between the loop integrals of the vector potential 
is a topologically invariant. We will show in this section 
that our discretized theory has the same property. 
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A. Commutators for the continuum case 


For the Chern-Simons gauge theory in the continuum, 
for two arbitrary curves C and C , we have the following 
commutation relation 



(7.1) 


where v[C, C'] is the number of (oriented) intersections 
between the two curves, i.e. the number of right-handed 
interactions of C and C' minus the number of left-handed 
onesP*^ 

If C and C' are closed loops, v[C, C] is a topologically 
invariant, and it is easy to verify that its value cannot 
change under any adiabatic procedures. In addition, if 
either C or C' can be contracted into a point (i.e. con¬ 
tractible), v[C, C'] = 0. 


B. Canonical quantization 


Using canonical quantization, it is straightforward to 
show that the conjugate field of the vector potential field 

A ei is 

= A K e>e ,A e , (7.2) 

oA e 27t 

This formula implies that for our discretized Chern- 
Simons theory, the vector potential A e (and linear su¬ 
perpositions of A e ’s) play both the role of the canoni¬ 
cal coordinates and that of the canonical momenta. Be¬ 
cause canonical coordinates and canonical momenta arise 
in pairs, this result requires that we must have even num¬ 
ber of linear independent A e s, i.e. the number of edges 
must be even. This is indeed true for any graphs con¬ 
sidered here. Utilizing the Euler characteristic, we know 
that the numbers of vertices, edges and faces must satisfy 
the following relation 


N v - N e + N f = 2 - 2 g, (7.3) 

where g is the genus of the underlying manifold. Because 
the vertex-face correspondence requires N v = Nf , the 
number of edge is 

N e = 2Nf - 2 + 2 g. (7.4) 

which is an even number. 

In canonical quantization, the commutator between 
a canonical coordinate and the corresponding canonical 
momentum is ih. Therefore, for our theory, we have 


Ae ’ 2n Ke '’ e " Ae 


= iS e 


(7.5) 


where 5 e y is the Kronecker delta and we set h to unity. 
Multiplying both sides by the inverse matrix of K , we 
obtain the commutation relation for the vector potential 


r 27T 1 27 T% i 

[A e , A e >] = i T K~} e = - —K~l, (7.6) 


Here, we used the fact that I\ 1 is an antisymmetric 
matrix. 

In order to ensure that the commutator [A e ,A e <] is 
nonsingular, we must require the K matrix being invert¬ 
ible. 


C. Paths, contractible and noncontractible cycles 


In this section, we will introduce two concepts from 
the graph theory: paths and cycles , which are discrete 
versions of curves and loops, respectively.^ 

A path is a sequence of vertices vq —> V\ —> Vi — > ... — > 
v m , in which any two consecutive vertices are connected 
by an edge. In the literature of graph theory it is often 
also assumed that a path never go through the same ver¬ 
tex twice. The length of a path is the total number of 
edges contained in the path. 

If vq v m , the path is called open. For vq = v m , the 
path is closed. A closed path (with nonzero length) is 
also called a cycle. In comparison with the continuum, 
it is easy to realize that open paths are discretized open 
curves, while cycles (i.e., closed paths) are discretized 
loops (i.e., closed curves). More precisely, a path (cycle) 
corresponds to a directed curve (loop), because a path 
(cycle) has a natural direction built in according to its 
definition, i.e. Vq v m . 

In the continuum, loops on a 2D manifold can be classi¬ 
fied into two categories: contractible or noncontractible, 
depending on whether or not the closed curve can be 
adiabatically contracted to a point. For a graph, there 
is a similar classification for cycles (closed paths) using 
a different but equivalent definition. We call a closed 
path (i.e. a cycle) contractible, if it is the boundary of 
some 2D area formed by a set of faces. Otherwise, it is 
noncontractible. For 2D closed and orientable surfaces 
in the continuum, noncontractible loops only exist for 
surfaces with nonzero genus (torus, double torus, etc.), 
while all loops on a genus zero surface (e.g. a sphere) 
are contractible. In graph theory, the same is true for 
cycles. For planar graphs defined on 2D closed and ori¬ 
entable surfaces, noncontractible cycle can only exist if 
the genus of the underlying 2D manifold is larger than 
zero. 

For a directed graph (or lattice), each path (P) can be 
represented by a W-dimensional vector, £p, whose eth 
component is 

! +l e £ P and e is along the direction of P 
— 1 e € P and e is opposite to the direction of P 
0 e^P 

(7.7) 

As will be shown in below, this object defines a dis¬ 
cretized line integral. In particular, if P is a cycle, 
provides a discretized loop integral. 
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D. Commutators and intersections 


For a path P on a graph G, we can define the integral 
(circulation) of the vector potential along this path as 


Wp = E,P,eA e 


(7.8) 


This object is the discretized version of a line integral 
f c A ■ dx along a path C. 

Now, we consider two different paths, P and P', and we 
define two integrals Wp and Wp< for P and P', respec¬ 
tively, using the definition of Eq. (7.8). In this section, 


we prove that the commutator between Wp and Wp< is 
determined by the number of oriented intersections be¬ 
tween the two paths v[P,P'\, 


C )'Tr r l 

[Wp,Wp-]=yP[P,P '] 1 


(7.9) 


which is the direct analog of the corresponding com¬ 
mutator of the Chern-Simons theory in the continuum, 
Eq. f7T). 


Utilizing the commutator of Eq. (7.6), we find 


27T2 

[Wp,Wp/] =-£p,e£p' ,e' 


-1 

e.e' 


(7.10) 


If the two paths P and P' share no common vertex, the 
intersection number is obviously zero v[P, P'] = 0. In the 
same time, \Wp,yV p>] als o vanishes, because every term 
on the r.h.s. of Eq. (7.10) is zero. 


If the two paths share some common vertices, only 
edges connected to thes e com mon vertices contribute to 
the commutator of Eq. 17. 101) , because K~ , = 0 for all 


other edges. Therefore, we only need to consider edges 
adjacent to each common vertex. As shown in Fig. [TJ 
we shall distinguish three different situations, shown in 
Figs. [Ml to |7 (c) | respectively, depending on whether 
the common vertex is a right-handed intersection, a left- 
handed intersection, or not an intersection. In Fig. [7] we 
label the edges of P as e\ and e 2 , whi le the edges of P' 
are called e[ and e 2 . Using Eq. (7.10), the commutator 
is given by 


[Wp,Wp>] ^ y: £p,ej£p%e'A'*. e / 

(=1 j=l 

= ~ ~Y2 EE^e i CP',e;')I ;ej ,e'»?2*;e i ,e' 

1=1 7=1 

(7.11) 

Here, we used the fact that K~ x = — K * and each el¬ 
ement of K* can be written as —r]lr] 2 /2 as defined in 
Sec. |VI[ For the first three figures in Fig. [TJ it is easy to 
verify that 


£p,e;£p',e'??2 ; ei) ej. 


+1 if i = j 
-1 if i ± j 


(7.12) 



(a)A right-handed intersection (b)A left-handed intersection 

(" = +!) (p = -U 



FIG. 7. (Color online) One common vertex shared by two 
paths. Here, we consider two paths P (thin red solid lines) 
and P' (thick blue solid lines). The arrows indicate the di¬ 
rection of each path. The disk in the middle is one common 
vertex shared by the two paths. Dashed lines represent other 
(possible) edges that are connected to the vertex, and they 
don’t contribte to the commutator that we want to compute. 
Figure (a) shows a right-handed intersection between P and 
P 1 and Fig. (b) is a left-handed one. In Fig. (c), the two 
paths don’t intersect. Figure (d) shows a special case, where 
one path terminates at this vertex. Here, the number of in¬ 
tersection can be ±1 or 0 depends on microscopic details. In 
Sec. |VII E[ a method will be introduced to obtain the value 
of v for Fig. (d) by defining a dual path in the dual graph. 


If the common vertex is a right-handed intersection 
of P and P', Fig. 7(a)| four possible cases need to 
be considered depending on the location of the face 
that paired up with the common vertex, i.e. (1) be¬ 
tween e! x and e 2 , (2) between e 2 and e' 2 , (3) between 
e' 2 and ei and (4) between ei and e' x . For case (1), we 

have = - 1 ’ dl, eiA = +1, <;e 2 ,e' = and 

r ll-e 2 ,e' 2 = -1 - For Case ( 2 )> we have = “l) 

^ ;ei ,e' = + 1 . ’tfiea.ei = + 1 and ^;e 2 l e^ = + L F ° r 
case (3), r i[. eu< = -1, vl, ei , e ' 2 = -1, ^ ;e2 ,e' = + 1 and 

^;e 2 ,e' = “ F For CaSe ( 4 )> = + 4 > ^ ;ei ,e' = + 1 ’ 

V*Ue 2 y, = + 1 and ^;e 2l e' = ~ F Usin S Ec l S - <0 


and (7.12), we find that for all these four cases, the com- 
mutator [Wp , Wp>] = 2ni/k. Therefore, we find that 
each right-handed intersection contribute 27 t i/k to the 
commutator. 

Using the same technique, we can prove [ yVp,Wp>] = 
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—27 ri/k for Fig. |7(b)| and [Wp,Wp'] = 0 for Fig. 7(c) 


In summary, we find that each right-handed (left-handed) 
intersection contribute +27 t i/k {—2iTi/k) to the commu¬ 
tator [Wp, Wp'], and thus we proved Eq. (7.9). 


E. Gauge invariance and the commutation relations 


In this section, we prove that the commutation rela¬ 
tions of Eq. (7.91 arise naturally, if we require the ac¬ 


tion to be gauge invariant, Eq. (4.5). In addition, a by¬ 


product of this proof offers a more rigorous definition for 
the number of oriented intersections, which eliminates 
the ambiguity demonstrated in Fig. |7(d)| There, the two 
path P and P' barely touch each other. Shall this counts 
as an intersection? This question will be answered in this 
section. 



encloses (and only encloses) these vertices (the dashed 
lines in Fig. [8]). This new cycle will be called the dual 
of P' and will be labeled as P'* . Here, we choose the 
direction of P'* such that its orientation is the same as 
that of P' . Below, we will prove that the gauge invariance 
immediately implies the commutator 

[Wp,WH = -^[P,P'*], (7.13) 


Here, instead of the number of intersections between P 
and P' , we shall count the number of intersections for 
P and P'*. Because the cycle P'* is defined in the dual 
graph, the number of intersections is always well-defined 
and this eliminates the ambiguity shown in Fig. |7(d)| 
Before proving Eq. (7.13), we would like to highlight 
that although the dual cycle P'* and the original cycle P 
are not identical, the difference between them is local and 
microscopic. This comes from the fact that our vertex- 
face correspondence is local , where a face is paired with 
one of its neighboring vertex. If we take the continuous 
limit and ignore differences at the microscopic level, the 
differences between P'* and P' vanishes, and therefore, 
we recover Eq. (7.1). 


Now we prove Eq. (7.13). First, we define a Nf- 
dimensional vector Qp ’ for the contractible cycle P', 
whose the /th component Qp>j is 


Qp’j 


1 if the face / is enclosed by P' 
0 if the face / is outside of P' 


(7.14) 


With this matrix Qp>j, the contractible cycle P' can be 
written as 


FIG. 8. (Color online) A cycle in a graph and the dual cycle 
in the dual graph. Here, we consider a planar graph with a 
local face-vertex correspondence. The vertices in the original 
graph is marked by disks, while the crosses label the faces, 
i.e. vertices of the dual graph. The local face-vertex corre¬ 
spondence is marked using the dotted lines, which pair up 
each face with one of its neighboring vertex. The thick solid 
(blue) lines marks a contractible cycle (a closed path) and 
the orientation of the cycle is marked by the arrows. For a 
contractible cycle, its interior is formed by a set of faces (dark 
region). For each face inside the dark region, we find the cor¬ 
responding vertex using the local face-vertex correspondence. 
These vertices are marked by circles. Then, we draw a loop 
in the dual graph, which encloses these vertices (red dashed 
lines connecting neighboring crosses). This loop in the dual 
graph is the dual of the original loop in the original graph. 
And we require the two loops to have the same orientation. 

Consider two paths P and P 7 . Here we assume that one 
of the paths is a contractible cycle (P 7 ), while the other 
is an open path with two open ends (P). As an example, 
a contractible cycle P' is plotted in Fig. [8] Because P 7 is 
contractible, it is the edge of an area formed by a set of 
faces (the dark region in Fig. [8]). Utilizing the vertex-face 
correspondence, this set of faces are mapped to a set of 
vertices, which are marked by circles in Fig. [8] Now we 
can define a cycle in the dual lattice such that the cycle 


£p',e — Qp>,f€f,e 


(7.15) 


wherejT e is defined in Eq. and £p ie is define in 

Eq. ( |7.7| ). The proof for Eq (7.15) is straightforward. 
For the r.h.s., it is easy to notice that for any e outside 
the region enclosed by the cycle P 7 , Qp',/£/, e = 0. For 
an edge inside the region enclosed by the cycle P 7 , it 
will induce two terms for the r.h.s., because each edge 
is shared by two faces. These two terms have opposite 
signs and thus cancel out, and thus <3p',/£/, e = 0. The 
only way to have a nonzero Qp',/£/, e is to require that 
e is an edge of the cycle P 7 . And it can be verified that 
the value and the sign of Qp',/£/, e match exactly £p', e - 
Similarly, in the dual space, we can write down the 
dual cycle P 7 * as, 


£Ke* = = -Q*P'., v D v , e 


(7.16) 


where Q* P „ /. = 1 for any faces (of the dual graph) inside 
the dual cycle P 7 *. Here, we relabeled the faces in the 
dual graph (/*) using corresponding vertices in the origi¬ 
nal graph ( v ). We also used the fact that = —D VtB 
as shown in Eq. (6.2). Because the vertices (of the orig¬ 


inal graph) enclosed by P 7 * are partners of the faces en¬ 
closed by P 7 , we have 


Q*p’* v = Qp\fM v j 


(7.17) 
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By combining the two equations above and relabeling e* 
as e, we find that 


£f>'*,e = v,e- 


(7.18) 


In the last step, we used the fact that the M matrix is 
orthogonal, M v j = MJ v . 


By substituting Eq. (7.15) into Eq. (7.10), we find that 


C )'TT1 

[Wp,W l P ]=-—^. e K-^ P ,. e 


= ~ -^£,P,e K e yQp'j£f,e' 

O.’JT't 

=- —£p ie Qp>jMjlD Vte . (7.19) 

Here, we utilized the condition of gauge invariance 
MylD Vte = Eq. ( |6.23| ). Using Eq. ( |7.18| ), the 

r.h.s. can be written as 


27TZ 


[Wp,W'p\ = -j-Sp,etp'-,e- 


(7.20) 


It is easy to verify that only intersections between P 
and P'* contribute to the r.h.s. of the equation. At a 
right-handed/left-handed intersection, £p ie £|>/» e = ±1, 
and thus 

\Wp,W p ] = ^v[P,P'*]. (7.21) 


VIII. WILSON LOOPS FOR 
NON-CONTRACTIBLE CYCLES 


We start this section by considering a planar graph 
embedded on a 2D torus (with genus g = 1). For this 
graph, there are two independent non-contractible cycles 
(i.e. discretized counterparts of the two non-contractible 
loops on a torus), which will be labeled as C and C' 
in this section. These two cycles intersect once with 
each other. Without loss of generality, we choose the 
oriented intersection number to be +1, instead of — 1, 


i.e., v[C,C'] = +1. As we proved above in Eq. (7.9), the 
commutator [Wc, Wc] = 2ni/k. 

Here we define Wilson loops for the two non- 
contractible cycles C and C' of the torus 


Wc = exp (iWc) 
Wc =exp (iW' c ) 


( 8 . 1 ) 

( 8 . 2 ) 


If we consider an eigenstate of Wc with eigenvalue w, 


W c |tf) = w |\k) 


(8.5) 


where w is a complex number, utilizing Eq. (8.4), it is 


straightforward to show that Wc I'F) is also an eigenstate 
of Wc and its eigenvalue is we~ 2m ' k 


W c (Wc\^)) = we- 2ni/k (Wc\^}) ( 8 . 6 ) 


In other words, we can consider Wc as a raising/lowering 
operator for the operator Wc, and vice versa. Starting 
from the eigenstate |\I/), eigenstates of Wc can be gener¬ 
ated by applying this raising/lowering operator, 

W c (Wg, I’E)) = we- 2nvi / k {W£, |T)) (8.7) 


i.e. Wq, 14') is an eigenstate with eigenvalue we ~ 2nm / k 

For an integer k , it is easy to note that when n = fc, the 
state Wq, | T) has the same eigenvalues as |\&). If |H>) 
and |T) are the same quantum state, Wp, I'll) generates k 
different eigenstates of Wp. From this results it follows 
the well known result that a Chern-Simons gauge the¬ 
ory has a fc-fold topological degeneracy on a torus. This 
conclusion is well known in the continuum. Our discus¬ 
sion above shows that the same is true in our discretized 
theory. 

It is straightforward to generate the discussion above 
to other 2D manifolds with different genus. For a pla¬ 
nar graph defined on a 2D surface with genus g , there 
are 2 g independent non-contractible cycles. As will be 
discussed in Sec. H (Fig. © , we can choose g of these 
cycles such that they don’t intersect with each other, Ci, 
C* 2 , ..., C g . The other g non-contractible cycles will be 
labeled C[, C ' 2 , ..., C' g . The absence of intersection for 
cycles Ct (i = 1,2,... ,g) implies that the Wilson loops 
defined on this cycles commute with each other, and thus 
we can consider common eigenstates for Wc x , Wc 2 , ■ ■ 
Wc g ■ The Wilson loops for the other g non-contractible 
cycles serve as raising and lowering operators. Starting 
from one common eigenstate of all Wc t s, we can use Wc 
to generate k 9 eigenstates (including the original one), 
which reflects the topological degeneracy of the Chern- 
Simons gauge theory on a surface with genus g. 


IX. LOCALITY 


Because the commutator [>Vc,W , c'] = 2ni/k is a com¬ 
plex number (i.e. is proportional to the identity opera¬ 
tor), it commutes with both Wc and W' c . Hence, using 
the Baker-Hausdorff-Campbell formula it follows that 

giVVcgiVVc _ giWcg!Wcg[!Wc,iW c /] (g g) 

and thus 

W c Wc = WcWce~ 2Ti/k (8.4) 


In this section, we verify that our theory is local. More 
precisely, we prove that (1) our action is local, i.e. the 
action does not have any coupling between fields that are 
not around the same face, (2) the flux attachment is local, 
i.e. for a charge at the vertex v. its magnetic flux must be 
located on a neighboring face, and (3) the commutator 
between vector potentials is local , i.e. for any two edges 
that do share a common vertex, the vector fields defined 
on them commute with each other. 
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A. The action 


In the discretized action of Eq. (3.4) there are no long 
range couplings beyond edg es an d vertices of the same 
face. In the first term in Eq. (3.4), because vanishes 
unless / and v are adjacent to each other, the action only 
contains couplings between nearby A v and A e (i.e. e and 
v must belong to the same face). For the second term, 
we know that K ete i = 0, if e and e! do not belong to the 
same face, and therefore, only short-range coupling (for 
edges of the same face) is included in this term. 


B. Flux attachment 


For the Chern-Simons gauge theory in the continuum, 
the flux attachment is local , i.e. for a point charge at Vq, 
the magnetic field is a delta function B oc 5(r — ro), and 
the B field vanishes when we more away from the point 
charge. 

In our discretized theory, this condition of locality is 
preserved to the maximum extent. As shown in Eq. (5.31, 
for a charge at the site v, the corresponding magnetic field 
only is present inside a single face , which is the closest 
analog of a delta function in a discrete setup. As for the 
relative locations of the charge and its flux, because these 
two objects on different parts of the graph (charges on 
vertices and fluxes on faces), it is impossible to require 
their location to coincide. Instead, we require the charge 
and the flux to be adjacent to each other. 

We emphasize that this locality condition plays a very 
important role, if we use the Chern-Simons gauge theory 
as a statistical field to change the statistics for matter 
fields coupled to it. To ensure that all particles have the 
correct statistics, when we move a particle A around an¬ 
other particle B , A must feel all the statistical field of B. 
In other words, no matter which path we choose, as long 
as A moves around B, the magnetic flux attached to B 
must be enclosed by the path of A. For our theory (and 
for the continuous Chern-Simons gauge theory), this is 
always true. However, if one were to violate the local¬ 
ity condition by putting the magnetic flux in a face not 
adjacent to the charge, it would be possible to move A 
around B without enclosing the flux inside the path. As 
a result, the statistics of the matter field would become 
ill defined. 


C. The commutation relations 


As shown in Eq. (7.6), the commutator for the vector 
potential is determined by the inverse of the K m atrix (or 
say the dual matrix K*). As proved in Sec. VI the K~ l 
matrix is also local, where the matrix element K~^, = 0, 


if e and e' does not share a common vertex. 

This results implies that a nonzero commutator can 
only arise for two neighboring edges, while for two edges 


separated away from each other (i.e. not sharing a com¬ 
mon vertex), the vector potential always commutes with 
each other. 


X. WHY N v = N f ? 


Above, we have shown that the existence of a local 
vertex-face correspondence is sufficient for the discretiza¬ 
tion of the Chern-Simons gauge theory. In this section, 
we prove that this condition is also necessary, if we want 
to preserve key properties of the Chern-Simons theory. 

Let us consider a generic discretized action of gauge 
fields A v and A e . Just as in the Chern-Simons gauge 
theory in the continuum, we assume that the action does 
not contain time derivatives of the time component of 
the gauge field A v , and that A v plays the role of a La¬ 
grange multiplier field that enforces a constraint on the 
local flux. For the coupling among the components of the 
gauge fields A e on different edges, the action only con¬ 
tains product between A e and d t A e f. We ignore possible 
terms with higher orders in time derivatives, which are 
less relevant in the sense of the renormalization group. 
In addition, we will only keep terms to the leading order 
in our action. 

With these assumptions, the most generic action that 
one can write down is 



dt 


A v M v ,f£ t f e A e ^A e K e yA e 


( 10 . 1 ) 


This action is very similar to the action we constructed 
above in Eq. (3.4). However, we must emphasize that 
here M and K are generic matrices, and that so far we are 
not putting any constraints on them. Most importantly, 
now we don’t require the graph to support a local vertex- 
face correspondence. Instead, we will consider generic 
situation and show that if we want the action to take 
this form, then the local vertex-face correspondence will 
arise naturally. 

In Sec. |X A[ we first introduce some mathematical tools 
from algebraic graph theory. Then, in Sec. |XB[ we will 
prove that the number of faces cannot exceed the num¬ 
ber of vertices (N v > Nf), and otherwise the theory will 
be singular. Then, in Sec. |X C[ we show that the flux 
attachment requires the number of vertices not to exceed 
the number of faces (N v < Nf). Combining these two 
conclusions together, we find that the graph must have 
the same number of vertices and faces (N v = Nf). Fi¬ 
nally, in Sec. |X D| we prove that a local vertex-flux cor¬ 
respondence is necessary, if we further require the flux 
attachment to be local. 


A. Edge-space, Cut-space and Loop-space 

Here we introduce some concepts from the algebra 
graph theoryj^that will be used later. In algebraic graph 
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theory, an A e -dimensional vector represents each edge e 
of a graph G, 


e e = (0,0,..., 1,..., 0), (10.2) 

where the eth component of the vector is 1 and all other 
components are 0. These vectors form the basis of a N e - 
dimensional linear space, which is called the edge space 
of the graph G. It is easy to realize that the A'-matrix 
defined above is a rank-2 tensor in this linear space. 

For a directed graph (i.e., a digraph), each (con¬ 
tractible or non-contractible) cycle C can be represented 
as an iV e -dimensional vector, £c, whose eth component 
is 

! +l e £ C and e is along the direction of C 
— 1 e £ C and e is opposite to the direction of C 
0 e^C 

(10.3) 


These vectors span a linear space, which is a subspace of 
the edge space. In algebraic graph theory, this subspace 
is known as the circuit-subspace. 

A cutset is a set of edges, where if we cut all the edges 
in a cutset, the graph is cut into two disconnected pieces. 
A more rigorous definition of a cutset relies on a partition 
of vertices. If V is the set of all vertices of a graph G, 
we can separate these vertices into two subset V-\ and V 2 , 
such that Vi U V -2 = V and V\ D V 2 = 0. This is called a 
partition of the set V. For each partition of V, we can 
define a cutset by collecting all edges of G that have one 
end in V\ and the other in V 2 . For a digraph, one can 
choose one of the two possible orientations for a cutset 
by specifying the vertices in Vi (or V 2 ) to be the positive 
ends, while the other to be negative. If an edge in the 
cutset points to the positive end of the cutset, it is a 
positive edge in this cutset. Otherwise, it is a negative 
edge. 

Similar to cycles discussed above, each cutset can also 
be represented by an W-dimensional vector £#, whose 
eth component is 

{ +1 e £ H and e is a positive edge 

— 1 e £ H and e is a negative edge (10.4) 

0 e^H 


The linear space spanned by these vectors is known as 
the cut-subspace, which is also a subspace of the edge 
space. For a planar graph, each cutset corresponds to a 
contractible cycle in the dual graph. 

In algebraic graph theory, it is shown that the edge 
space is the direct sum of the circuit-subspace and the 
cut-subspace. In Appendix [F] we provide a proof for the 
planar graphs considered here. This result implies that 
for the edge space, instead of using the basis shown above 
in Eq. (10.2), we can choose a new basis for the edge space 
by selecting a complete basis of the circuit-subspace and 
a complete basis of the cut-subspace. 

For planar graphs, we can use all independent (con¬ 
tractible or noncontractible) cycles to form a basis for the 


circuit-subspace. For the cut-subspace, all independent 
contractible cycles in the dual graph forms a complete 
basis. Therefore, we can span the edge space using these 
loops. Using this new basis, we can rewrite all tensors de¬ 
fined on the edge space, including the A -1 matrix, which 
will be done in the next section. 


B. N v > N f 


We will now prove that for the A matrix to be nonsin¬ 
gular and the discretized t heory to preserve the correct 
commutation relation of Eq (7.1), the number of faces can 
never exceed the number of vertices. Using the generic 
action shown in Eq. (10.1) (remember that A and M are 


now two arbitrary matrices), we find that for the generic 
setup, we shall still expect the commutation relation 


[A e ,A e ,] = ~K~l, 


(10.5) 


Because singularities in the commutation relations must 
be avoided, the K matrix must be invertible. In addition, 
if we consider two cycles (loops) G and C', we shall expect 
the commutation relation 

O7 rn 

[W c , Wc'] = —j£~ v [G, C'] , (10.6) 

As shown above, this commutator is a topological invari¬ 
ant and it is one of the key feature of the Chern-Simons 
gauge theory. Thus, we will require Eq. (10.6 1 for our 
discretized theory. 

Below, we prove that if we assume the topologically 
invariant commutation relation, Eq. (10.6), then the K 


matrix must be singular if N v < Nf. Therefore, we must 
have N v > Nf. We will start from a genus zero surface 
and then expand the conclusion to other surfaces with 
nonzero genus. 


1. Graphs on a genus zero surface 


Here, we consider graphs defined on a genus zero sur¬ 
face (a sphere). Instead of directly showing that the K 
matrix is singular for N v < Nf, here we take a different 
but equivalent approach. We will start by assuming the 
K matrix is invertible and work with the A -1 matrix. 
Then, using the commutation relation, we will show that 
the determinant of A' -1 matrix is zero for N v < Nf, and 
thus the K matrix is singular. 

Using Eqs. (10.5) and (10.6), we know that 


K e i£,C,e£,C',e 


= -u[C,C'} 


(10.7) 


Here, we choose a new basis set for the e dge s pace. 
Instead of using the vectors shown in Eq. (10.2), we 
use a set of vectors & with i = 1,2 ,...,N e . For 
i = 1, 2,... Nf — 1, are independent cycles, i.e. they 
form a complete basis of the circuit-subspace. For i = 
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FIG. 9. (Color online) Non-contractible cycles on a surface 
with nonzero genus. For a genus g surface, we can choose g 
independent non-contractible cycles, which do not intersect 
with one another. These g cycles will be used as vector 
with i = Nf,Nf + 1,..., Nf + g — 1, in our complete basis 
for the edge space. Here, we show an example with g = 3. 
The three red loops marks three independent noncontractible 
cycles without any intersections. 


Nf,Nf + l,...,N e , the corresponding 6 are independent 
cutsets, i.e. they are a complete basis of the cut-subspace. 
Using this new basis, we can define a K -1 matrix as 


—l 


Kr / = K p . 

L.j 


'Ci,e£j,e 


( 10 . 8 ) 


For i and j smaller than Nf, and £j are contractible 
cycles of the graph (for a planar defined on a closed ori- 
entable 2D surface with genus zero, all cycles are con¬ 
tractible). Using Eq. (10.7), it is easy to realize that 
K~j = 0 for i. and j smaller than Nf. (As shown above, 
the number of oriented intersection for contractible loops 
is always zero). Therefore, we can write the A " -1 matrix 
in a block form 


= (-^ b) < ia9 > 

Here the first block 0 is a (Nf — 1) x (Nf — 1) zero matrix 
and B is a (N e — Nf + 1) x (N e — Nf + 1) matrix. Using the 
Euler characteristic N v — N e +Nf = 2 — 2 g, we can rewrite 
the dimensions of B as (N v — 1) x (N v — 1), since we have 
assumed the genus being zero, g = 0. The block A has 
dimension(7Vy — 1) x (N v — 1) and A T is the transpose of 
A. 

For a matrix with a block of zeros as shown in 
Eq. ( |10.9| ), the determinant of the matrix must be zero, 
if the zero block is larger than the B block (see Appendix 
Sec. [g] for a proof). Therefore, if N v < Nf, det A'” 1 = 0. 
Because £, is a complete basis for the edge space, this im¬ 
plies that det A" -1 = 0 and thus A' is a singular matrix. 


2. Surfaces with nonzero genus 

For a surface with nonzero genus, the same conclusion 
can be proved. Here, we choose the following basis of 
the edge space £ 1 , 62 , •••&v e - For i = 1 , 2 ,... IV/ — 1 , 
6 are Nf — 1 independent contractible cycles. Then, 
for i = Nf, Nf + 1,..., Nf + 2g — 1, 6 are independent 
noncontractible cycles. For these noncontractible cycles, 


we choose to have the first g noncontractible cycles (Nf < 
i < Nf +g — 1) not to cross with one another as shown in 
Fig. [9} It is easy to realize that the first Nf+2g— 1 vectors 
here form a basis of the circuit subspace, while the rest 
are chosen to be a complete basis of the cut-subspace. 

Using this new basis, we can define a AT -1 matrix as 


K^ 1 = A '” 1 £■ £ / 


( 10 . 10 ) 


For i < Nf and j < Nf + g (or i < Nf + g and i < Nf), 
K~J = 0. This is because here 6 and £.,■ are to cycles 
of the graph G and at least one of them is contractible. 
According to Eq. (10.7), A'T 1 = 0, because the number 
of oriented intersections vanish when one of the cycle is 
contractible. For Nf < i < Nf + g — 1 and Nf < j < 
Nf + g — 1, 6 and £j are two noncontractible cycles, 
but we have required that these cycles do not cross each 
other, i.e., zz[£j,£j] = 0, and thus, K~j = 0. 

With this knowledge, we can write the matrix Kj j in 
this block form 


Kij = f) ( 10 . 11 ) 


Here, the upper-left Conner to be a zero matrix with di¬ 
mension (Nf+g— l)x (Nf+g— 1). The dimension of the B 
matrix is (N e — Nf—g+l)x(N e —Nf—g+l). Utilizing the 
Euler characteristic N v — N e + Nf = 2 — 2 g, we find that 
the dimensions of B is in fact (N v + g — 1) x (N v + g — 1). 

If Nf > N v , again, we find that the 0 block is larger 
than the block of B, and therefore, det A " 1 = 0 (see Ap¬ 
pendix [G] for a proof). Because { 6 } is a compete basis 
for the edge space, this implies that det A '” 1 = 0, and 
thus AT is a singular matrix. 

In summary, we proved that in order to preserve the 
commutation relations, Eq. (10.6), we must have N v > 
Nf. Otherwise the K matrix would be singular. 


C. Flux attachment and N v <N f 


Let us now prove that the flux attachment also requires 
N v < Nf. Flux attachment implies that for each charge 
distribution, there is a corresponding unique distribution 
for magnetic fluxes. Because charge can be distributed 
on N v sites, to ensure that there is a corresponding flux 
distribution for every charge configuration, we must have 
equal number or more faces to put the fluxes. 

A more rigorous proof can be formulated by taking 
a functional derivative to the generic action Eq. (10.1), 


SS/SA V , which result in the flux attachment condition 


Q v 2n M vJ * f , 


( 10 . 12 ) 


If we want the flux attached to a charge to be local (i.e. 
the flux for a point charge occupies only a single face), 
for each vertex v, the M v j is nonzero only for one value 
of /. As a result, the M matrix defines a mapping from 
v to /. 
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This mapping must be injective. Namely, for two dif¬ 
ferent vertices, their corresponding faces must be differ¬ 
ent. This is so because if two different vertices v and 
v' are mapped to the same face /, then Eq. ( 10.12[ ) will 
require that q v = q ' v , i.e. two different vertices always 
have the same charge, which is obviously not a physi¬ 
cally necessary constraint. Thus, for a injective mapping 
from vertices to faces, we must have N v < Nf. 


tetrahedron is self dual (i.e. the dual graph is also a tetra¬ 
hedron). This is also the simplest setup for discretizing 
the Chern-Simons gauge theory. 


A. The action 


D. local vertex-surface correspondence 


We label the vertices, edges and faces of a tetrahedron 
as shown in Fig. |10| In this convention, the incident 
matrix D v>e [Eq. (|4.3|] is 


In the previous two subsections we proved that N v < 
Nf and N v > Nf must hold simultaneously. Therefore, 
the graph must have the same numbers of vertices and 
faces N v = Nf. With N v = Nf, the mapping from ver¬ 
tices to faces discussed above become a one-to-one cor¬ 
respondence between vertices and faces. As addressed in 
Sec. |IX[ it is important to ensure that this correspon¬ 
dence is local. As a result, the local vertex-face corre¬ 
spondence arises naturally, when we try to ensure the 
theory being nonsingular and the key properties of the 
Chern-Simons gauge theory is preserved. 


/ 0 +1 -1 +1 0 0 \ 
- 10 + 10+101 
+ 1-10 0 0+1 
Vo 0 0-1-1 -1/ 


and the £/ je matrix defined in Eq. (3.6) is 


/+1 0 0 0 +1 - 1 \ 
0+10-10+11 
0 0+1+1-10 
V-l -1 -1 0 0 0 / 


( 11 . 1 ) 


( 11 . 2 ) 


XI. THE CHERN-SIMONS GAUGE THEORY 
ON A TETRAHEDRON 



FIG. 10. (Color online) A tetrahedron viewing from the top. 
The circles represent the vertices, and the arrows are the edges 
(direction assigned). Here, we label each vertices, edges and 
faces using integers. The last face (i.e. the face 4) is on the 
other side of the tetrahedron invisible from the current view 
point. 


In this section, we demonstrate our generic theory by 
presenting a specific example, i.e. by discretizing the 
Chern-Simons gauge theory on a tetrahedron. A tetra¬ 
hedron is a planar graph defined on a manifold with g = 0 
(a sphere). In addition, it is easy to verify that a tetra¬ 
hedron satisfies the cretieron presented in Sec. [TTJ and 
thus a discretized Chern-Simons gauge theory can be 
constructed. It is also worthwhile to emphasize that a 


We choose the local vertex-face correspondence such 
that vertices 1, 2, 3 and 4 pairs up with faces 2, 3, 4 and 
1 respectively. Therefore, the M v j matrix is 


M = 


/0 1 0 0\ 
0 0 101 
0 0 0 1 
Vl 0 0 0/ 


(11.3) 


Using this vertex-face correspondence, we can get 
the K e y matrix following the procedure described in 
Sec. |III B[ which is 


/O 

+i 

+1 

0 

-1 

- 1 \ 


-1 

0 

-1 

+1 

0 

-1 

1 

-1 

+i 

0 

-1 

+1 

0 

2 

0 

-l 

+1 

0 

+1 

-1 


+1 

0 

-1 

-1 

0 

-1 


V+i 

+1 

0 

+1 

+1 

0 / 


Using these two matrices, the action can be written down 
as shown in Eq. (3.4). 


Here, we can verify easily that the I\ matrix is in¬ 
vertible. In addition, it is straightforward to show that 
M t = M~ x and K T = —K (i.e. M is an orthogonal 
matrix and K is antisymmetric), in agreement with the 
generic result proved above. 


In addition, it is also straightforward to verify th at th e 
matrices satisfy the gauge invariance condition, Eq. (4.51, 
because = DK T . 
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B. Dual graph 


In the dual graph, it is straightforward to get the dual 
of the incident matrix and that of the £ matrix. 


D* = 


+ 1-10 0 0 +1 — T \ 

0+10-10+1 
0 0+1+1-10 
V-1 -1-10 0 0 


(11.5) 


and 


r = 


0 -1 +1 -1 0 

+ 10-10 -1 
-1+10 0 0-1 
0 0 0 +1 +1 + 1 / 


( 11 . 6 ) 


By comparing the D and £* (£ and D*) matrices, we 
find that D* = £ and £* = —ID, which verify Eqs. (6.1) 
and (6.2). 


In the dual graph, if we use the same vertex-face cor¬ 
respondence, we get the M* matrix 


M* = 


/0 0 0 h 
10 0 0 
0 10 0 
Vo 0 1 0/ 


(11.7) 


and the K* matrix is 


K* 



(° 

+1 

+1 

0 

-1 

- 1 \ 


-1 

0 

-1 

+1 

0 

-1 

1 

-1 

+1 

0 

-1 

+1 

0 

2 

0 

-1 

+1 

0 

+1 

-1 


+1 

0 

-1 

-1 

0 

-1 


\+l 

+1 

0 

+1 

+1 

0 / 


( 11 . 8 ) 


Using the matrices M*, K* and £*, we can write down 
the action for the discretized Chern-Simons gauge theory 
in the dual graph using Eq. (6.3). 

Here, we can verify that K* = —A' -1 and M* = 
M -1 , as well as the gauge invariance condition M*£* = 
D*(K *) t , Eq. (pUnl 


XII. CONCLUSIONS AND DISCUSSION 

In this paper, we proved that the Chern-Simons gauge 
theory can be discretized for generic planar graphs on 
arbitrary 2D closed orientable manifold as long as a local 
vertex-face correspondence can be defined on the graph. 
This condition is also necessary, if we want the theory 
to be nonsingular and to preserve some key properties of 
the Chern-Simons gauge theory. In particular, we showed 
that the gauge invariance of the discretized theory re¬ 
quires that the vertex-face correspondence to be strictly 
enforced. 

We also find a necessary and sufficient condition, which 
an be used to determine whether such a correspondence 
can be defined on a particle graph or not, based on the 


number of faces and vertices in this graph and its sub¬ 
graphs. 

The generalized discretized Chern-Simons gauge the¬ 
ory that we presented here has a number of interest¬ 
ing applications. One direction of further research is 
to consider the fractional quantum Hall effect on lat¬ 
tices, a problem that has not attracted much attention so 
f ar . HUM A more general theory of the fractional quantum 
Hall effect on lattices is of interest in the context of frac¬ 
tionalized time-reversal breaking top ological insulators so 
far as adiabatic continuity holds“*2“ These methods are 
also relevant to frustrated quantum antiferromagnets as 
we showed recently.^ 

There are several open as yet unsolved issues. One is 
to relax somewhat the vertex-face correspondence. Since, 
as we showed, this is required by gauge invariance, any 
violation of this correspondence is equivalent to either 
the insertion of background static charges or background 
static fluxes. This viewpoint may offer a way to gener¬ 
alize this construction to other lattices (e.g. triangular 
and honeycomb) as well as to investigate the role of lat¬ 
tice topological defects such as dislocations and disclina- 
tions of time-reversal breaking fluids, including quantum 
Hall fluids, where the role of geometry has been focus of 
recent interest EDED 

As a side comment, it is also worthwhile to note that 
two of the graphs shown in Fig. [2] (the Kagome lattice 
and the dice lattice) belong to the family of isostatic 
lattices. The terminology of isostatic lattices is devel¬ 
oped in the study of mechanical stability transition^ 
and recently, topologically nontrivial elastic modes are 
observed in some of these isostatic systems, including 
protected zero-energy edges states, nontrivial topologi¬ 
cal indices and topological zero-energy solitonJ^MH] Al¬ 
though the topological nature of those isostatic elastic 
systems are very different from a Chern-Simons gauge 
theory, it is not an accident that same lattices arises in 
these two seemingly unrelated areas. As shown in Ap¬ 
pendix (Sec. the isostatic condition is closely related 
with (and slightly stronger than) the criterion for the ex¬ 
istence of local vertex-face correspondence, which is the 
fundamental reason why some lattices can be used for 
both studies. 
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Appendix A: Simple graph 


Here, we demonstrate the definition of a simple graph 
by presenting situations that are not allowed in a simple 
graph, as shown in Fig. ED 



FIG. 11. Examples of graphs that are not simple. The fig¬ 
ures demonstrate situations that are not allowed for a simple 
graph. Figure (a), shows a pair of sites connected by three 
different edges. In Fig. (b), one of the edges connects a site 
with itself (i.e. the two ends of a edge coincide). 


Appendix B: Local vertex-face correspondence 

In this section,we prove that criterion presented in 
Sec. [TT] is a sufficient and necessary condition for a graph 
to have a local vertex-face correspondence by mapping 
this problem to Hall’s marriage problem 47 . 

The marriage problem considers a finite set of girls, 
each of whom knows several boys, and the task is to find 
the sufficient and necessary condition, under which all 
the girls can marry the boys in such a way that each 
girl marries a boy she knows (marriage is assumed to be 
one-to-one here). The solution to the marriage problem 
lies in Hall’s marriage theorem , which states that a nec¬ 
essary and sufficient condition for the existence of such a 
matching is that each set of k girls collectively knows at 
least k boys. 

Our goal here is to identify the sufficient and necessary 
for a graph satisfies the local-flux-attachment condition. 
And this problem can be mapped to the marriage prob¬ 
lem by mapping faces into girls and vertices as boys. If 
the vertex v is the vertex of the face /, the corresponding 
boy and girl know each other. Under this mapping, the 
local-flux-attachment condition is exactly the marriage 


problem. 

Here, we first prove that the criterion presented in 
Sec. [II] is a necessary condition. We consider a subgraph 
of G. By mapping to the marriage problem, all the faces 
in the subgraph forms a subset of girls, and the boys that 
they know are included by the set of vertices of the sub¬ 
graph, i.e. if the subgraph contains Nj faces (girls) and 
TV' vertices (boys), the number of boys that these girls 
know is equal to or smaller than N' v . Therefore, based 
on Hall’s marriage theorem, we must have TV' > N'j in 
every subgraph, if every face can marry a vertex that is 
adjacent to it. In other words, this criterion is necessary 
for the existence of a local vertex-face correspondence. 

We can also prove that the criterion is sufficient by 
considering subgraphs that satisfying the following con¬ 
dition: every vertex in the subgraph is adjacent to at 
least one face of the subgraph. (This condition does not 
hold for all subgraphs. For example, if a subgraph con¬ 
tains dangling bonds, the vertex located at the free end 
of a dangling bond is not adjacent to any faces in the 
subgraph). For these subgraphs (in which all vertices are 
adjacent to at least one face of the subgraph), the num¬ 
ber of vertices (IV') equals to the number of boys that are 
known by the girls (faces) in the subgraph. And there¬ 
fore, if the criterion in Sec. [TT] is satisfied, the marriage 
theorem ensures immediately the existence of (at least) 
one local vertex-face correspondence. 


Appendix C: Gauge symmetry 


Here, we prove that Eq. (4.5) is the sufficient and nec¬ 
essary condition to maintain the gauge symmetry in our 
theory [Eq. fl3l] ) ] 

Fi rst, we substitute the magnetic flux in Eq. (3.4) by 
Eq. pfil), 


S — —— [ dt(A v M vJ £ fte A e — -A ei K e e >A e f') (Cl) 


2-7T 


Under the gauge transformation 


A v —>A V — do (f> v 
A e >A e D ve cf) v 


the action Eq (Cl) is transfered to 


(C2) 

(C3) 


S^S + 


2tt 


dt 


(pvMv,f£,f,eAe T ^D V}f ,(j) v K e ^ e A e i -f- ^ A e K e e tD v e r 


) J dt ( Dv ’ e< t >vKe ’ e ' Dv '’ e "i )v '') ( C4 ) 


Here, the second term on the r.h.s. is linear in <f>, while 
the last term is 0((j) 2 ). In order to preserve the gauge 


symmetry, we need both these two terms to vanish, i.e. 



Mvj£f,eA e — A e K e ^D Vfi 


dt 





= 0 
= 0 


(C5) 

(C6) 
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In Eq. (C51, we used the fact that 


dtD v ^ e (j) v K e ^A e i I dtA e K eie iD Vie f(j) v 


(C7) 


which can be proved via integrating by part and realizing 
that the K matrix is antisymmetric. 

Equations (C5) and (C 61 imply that 


(C 8 ) 

(C9) 


— -^e,e' D v ,e' 

D v ' ,e' = 0 

These two conditions are not independent to each other. 
In fact, Eq. (C 8 ) automatically implies Eq. (C9). Thi s 
can be realized by noticing that according to Eq. (C 8 ), 
we have 


Dy_ oKf>_pf D v 




(CIO) 


The r.h.s. of this equation is zero because D v e ^f^ e = 0, 
and thus Eq. (C91 arises automatically. 


Here, we explain why e £/, e = 0. For any fixed /, 
£/ je represent an loop in the graph. If v is not a vertex on 
this loop, D Vte t;f' e must vanish, because D v e = 0. If the 
loop paths through v , there must be two edges along these 
loop that are connected to v. which we will call e\ and 
It is easy to realize that according to the definition of £ 
and D, D„ iei £/ iei = —D Vt e 2 ^f^ 2 (here, we don’t sum over 
repeated indices e\ and e^)- Therefore, the contributions 
to D v e ^f e cancels out, i.e. D v ^f, e = 0. This relation 
can also be written in a matrix form and the same is true 
for the dual graph 


D( r = D*(C) T = 0 (Cll) 

which will be used below in Sec. [e] Here £ T represents 
the transpose matrix of £. 


Appendix D: The directions of edges 


In this sectio n, w e prove that the condition of gauge 
invariance [Eq. (4.5)] is independent of the choice on the 
edge directions. 

As shown in the main text, we assign a direction for 
each edge in order to define the vector potential on a 
graph. These directions can be assigned in arbitrary ways 
and the choice of directions will not have any impact 
for any physics properties. For the condition of gauge 
invariance [Eq. (4.5)], this statement is also true. 

To prove this statement, we flip the direction of an 
arbitrary edge eo and consider two different situations 
eo = e and eo A e. 

If e 0 = e, as we fli p th e direction assigned to the edge 
eo, the l.h.s. of Eq. ( |4.5| ) changes sign, because M v j —>• 
M v j and £/ j(3o —> — £/, eo . The r.h.s. of the equation 
also flips sign, since K eo y —> —K e0tB > and D v e > —> D ve <. 
Because both sides of the equation flips sign when we 
flip the direction of eo, the equation remains invariant 


and thus is independent of the choice of the direction of 
eo 

If e 7 ^ eo, the l.h.s. of Eq. (4.5) remains invariant, 
because neither M v j nor £j je relies on the direction of 
eo- For the r.h.s., because both A' e>eo and D eo>v flip signs 
as we flip the direction of eo, their product remains the 
same. As a result, the equation is again independent of 
the direction of eo- 

This conclusion implies that in order to prove Eq. (3.6 1 , 
it is sufficient to verify the formula for just one specific 
choice of edge directions. 


Appendix E: Lattice duality 


In this section, we prove that the discretized Chern- 
Sirnons theory on the original lattice and the dual lat¬ 
tice are dual to each other by coupling the Chern-Simons 
gauge theory with gauge fields on the dual lattice. 


S = Scs + ^coupling 


(El) 


Here, the first term is our discrete Chern-Simons gauge 
theory 

Scs = 2 - J (^A v M v j$f — -A e K e yA e ^j (E 2 ) 

and the second term couples the Chern-Simons field A 
with gauge fields on the dual graph a* 

-^coupling = [ ^{Cf,e* a *e* A vSf‘,v + e ,a* v ,A e 5 e ^ 


- d 0 a* e ,A e 5 ete *) 


(E3) 


As will be shown in Sec. |E 1 [ this coupling is gauge in¬ 
variant. 


Below, in Sec. E 1 we first prove that same as in the 


continuum, the dual gauge field a* can be used to de¬ 
scribe the charge and current on the original lattice and 
we will also show that 5 coup ii ng is gauge invariant. Then, 
in Sec. |E 2[ we show that by integrating out the A field, 
the dual theory is obtained, which matches exactly the 
discrete Chern-Simons field on the dual lattice (but with 
a different coupling constant k* = — 1 /k). Because our 
action describes a quadratic theory, this calculation is 
exact. 


1. Gauge field on the dual lattice 

Same as in the continuum, we can consider the dual 
gauge field a* (defined on the dual lattice) as a descrip¬ 
tion for the charge and current on the original lattice. 
Here, the charge that resides at each vertex is called p v 
and the current on each edge is labeled as j e . 


1 /-* * 

Pv “ 2 n^ f *’ e * ae * 

(E4) 

j e (D*, e , a *, - d 0 a*,) 

(E5) 
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In Eq. (E4), we choose f* = v, and for Eq. (E5|, e = e*. 


It is easy to verify that the charge and current are gauge 
invariant and satisfy the charge conservation law (i.e. the 
continuity equation) 


&tPv k) v ,eje d 


(E6) 


where dt is the time derivative and D v e (the incident ma¬ 
trix) plays the role of a discretized divergence (multiplied 
by -1). 

To prove the gauge invariance, we perform the gauge 
transformation for a* 


a*. —>a*» — do <t>*. 

a* —>-a*. - D* .. 


Under this transformation, p v and j e are invariant 


pv ~*Pv — £f*,e*D. 


v* ,e* tv* 


— Pv 


(E7) 

(E8) 


(E9) 


Je tje - ^ {D* v . <e .dof v . - doDl.^l,) = j e (E10) 


For Eq. (|E9j), we used the fact e »U*» e , = 0. This 
relation is proved in Eq. (Cll) and it i s a d iscretized ver¬ 
sion of the formula VxV = 0. In Eq. (|E10|), 


and —doD*, e *<j>** cancels out because the incident ma¬ 
trix D*, e * is time-independent and thus commute with 
the time derivative <9 0 - 

The continuity equation can be proved using the fol¬ 
lowing two equations 


do Pv = a*. 


(Ell) 


2n ' 


(E12) 


In the second equation here, we used the fact that D v e = 
[Eq- f6|l)] and ^. >e .= 0 [Eq. ((ClT])].’ By 


adding the two questions together, the continuity equa¬ 
tion is obtained. 


Using Eqs. (E4| and (E5), the coupling between the A 
and a* fields [i.e. Eq. (E3iJ can be rewritten as 


^coupling — J dt (p v A v -f- j e A e ) . (E13) 

Because both p v and j e are gauge invariant, the coupling 
term must also be gauge invariant. 

The coupling S^upiing is also invariant under gauge 
transformation 


A v ~^A V — do4>v 
A e y A e D v e <j> v 


(E14) 

(E15) 


Using Eq. (E13), we find that the gauge transformation 


-A,,eJe = if},, e ,- doal.) = -~e f . ie .doa* e . 


turns ^coupling into 

^coupling — ^coupling j dt (p v (){)(p v T jeD l; e ) 

= ^coupling T J dt (doPv jedd v ,e) 4*v (E16) 

After a integration by part (for t), the last term in this 
formula vanishes due to the continuity equation, and thus 
^coupling is gauge invariant. 

2. Duality transformation 

In the path integral approach, a gauge fixing term 
•Sgauge fixing needs to be introduced 

*~*CS 4“ ^coupling 4“ ^gauge fixing 

Without loss of generality, here we choose 


^gauge fixing 


[ dt 1 

( dA v 

dA v \ 

/ 2tt ' 

\ dt 

dt J 


(E17) 


(E18) 


In the frequency space, the action of Eq. (E171 takes the 
following form 


S =Scs 4- ^coupling *$gauge fixing 


= A e (w)) 

OJ 


( auj 2 /k 

\m) T 


Mi \ 
—iioK J 





( A v (-w)\ 
\A e (-u))J 


(E19) 

(E20) 


Here, we write the Lagrangian as block matrices. Bold letters in this equation are vectors. For example, a*» represents 
a N e * -dimensional vector, whose components are a*, on each edge. The same is true for a*„, A e , or A v . The first 
matrix in the equation above contains Scs and Sgauge fixing, while the second matrix is for Sampling- 

By integrating out the A field, we obtain a dual gauge theory for a* on the dual graph. For a quadratic theory as 
shown above, this can be done exactly. 


S 


1 

2k 


«* ( w ) 


a e* M) 


GcV 



(atu 2 /k Mi \ 1 / 0 

{(Mi) T -iuK) -zcoj 


(E21) 
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The inverse matrix in the equation above can be computed using the technique of blockwise inversion 


A B 
C D 


{A- BD^C)- 1 -(A - BD~ 1 C)~ 1 BD~ 1 

-D~ l C{A - BD-'C)- 1 D- 1 + D~ 1 C{A - BD~ 1 C)~ 1 BD~ 1 


(E22) 


where A, B , C and D are matrix sub-blocks. For our matrix inverse, the block A is an identity matrix (multiply 
by a real number aw 2 /fc), and we have B = C T = and the block D is —iivK. Using the commutation relation 
Eq (7.10) and the fact that this commutator is zero for two contractible loops, it is easy to show that 


=o 


(E23) 


i.e., BD l C = 0. Therefore, we find that 

-l 


(au 2 /k 


\(M£) T —icoK) 


— %K~ i eM T -Aif - 1 - K~ 1 £ t M t M£K - 1 


— %M£K~ 

°j^ J ’ 


Now, we will use Eq. (6.20), which tells that 


M^IC- 1 =-D = £* 
k~H t m t = D t = -(D t 


(E24) 


(E25) 

(E26) 


Here, we also used the fact that K is an anti-symmetric matrix and Eq. (6.2) (D = —£*). Using these two relations, 
we find that 


(aco 2 /k 

l (M£) t -iojK 


-l 


fc _ ik c* 

otto 2 czuj 3 y 

\(C) T -±K- 1 + Ar(D T r 


(E27) 


And therefore, 
1 


5 _ ~ E ( a v* H a e* M) 


ICJ 


ik e* 
aw 3 ^ 


lk (n T + 

a;,(-w) 
(■ D*) T -iuj U£.(-w) 


r \ ( a*»(—w) 


= -^EK-o 


w a_, (w 


)) 


0 -±D*K- X 

IUJ 

0 -AT" 1 




0 


* ir-l 


£>*AT 


-K~ 1 (D*) J iuiK 


'—i 


a e* (~W) 


(E28) 


Here, we used the fact that D*(£*) T = 0 [Eq. (Cll)] and D*I\ 1 {D*) T = £K 1 £ r = 0 [Eqs. (6.1) and (E23)]. 
Because K~ l = —K* and D*K* = — M*£*, we get 


4E(<*H a e*M)L 


2k 


(D*K*) t —iujK* J la *.(-w) 


-me 


a*,(-u) 


= - s EM«) *•(")) 


0 M*r \ /a*,(- W ) 
(M*C) t -iujK* la*.(-w) 


(E29) 


By transferring from the frequency space w back to time 
t, we find 


5 = 


-1/fc 

27T 




<,M*. - —a*mK*» t 


(E30) 


This is exactly our discrete Chern-Simons gauge theory 
defined on the dual graph [Eq. (6.3)] with topological 
index k* = —1/fc. 


Appendix F: Edge space, circuit-subspace and the 
cut-subspace 


In this section, we prove that the edge space is the di¬ 
rect sum of the circuit-subspace and the cut-subspace. 
Although this conclusion applies generically to planar 
and non-planar graphs, we will only discuss planar graphs 
here for simplicity, since the manuscript only consider 
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planar ones. 

We first prove that the circuit-subspace and the cut- 
subspace are orthogonal to each other. This can be veri¬ 
fied easily by noticing that vectors from these two spaces 
(£c and £h) are orthogonal to each other (i.e. their dot 
product is zero) 

£c,e,£ir,e = o. (FI) 

In addition, we can prove that the dimension of the 
circuit-subspace plus the dimension of the cut-subspace 
coincides with the dimension of the edge space. Com¬ 
bined with the orthogonality proved above, this conclu¬ 
sion implies that the direct sum of the circuit-subspace 
and the cut-subspace is the edge space. 

As mentioned in the main text, the basis of the circuit- 
subspace can be formed by all independent (contractible 
and non-contractible) cycles of a graph. For a planar 
graph with Nf faces defined on a manifold with genus 
g , there are Nf — 1 independent contractible cycles and 
2 g independent non-contractible cycles, i.e. Nf — 1 + 2g 
independent loops in total. Therefore, the dimensionality 
of the of circuit-subspace is Nf — 1 + 2 g. 

For a planar graph G, cutsets corresponds to con¬ 
tractible cycles in the dual graph G*. Because the dual 
graph has N v faces, same as the number of vertices in 
the original graph, the number of independent cutsets 
(i.e. the number of independent contractible cycles in 
the dual graph) is N v — 1. 

If we add the dimensions of the cut subspace and the 
circuit subspace together, we get Nf+N v — 2 + 2g, which 
coincides with the number of edges N e . i.e. the dimension 
of the edge space. Here we utilized the fact that a closed 
orientable surface with genus g , the Euler characteristic 
is 2 — 2<7 

N v - N e + Nf = 2 - 2 g. (F2) 

Because the circuit-subspace and the cut-subspace are 
two orthogonal subspaces of the edge space, and the total 
dimensions of these two subspaces match the dimensions 
of the edge space, we proved that the direct sum of these 
two subspace is the edge space. 

Appendix G: the Determinant of a block matrix 

Here, we consider a ( N + M) x (N + M) matrix 

M = (d b) < G1 > 

where each letter in the matrix represents a block matrix 
and the 0 ( B ) matrix has dimensions N x N (M x M). 
We will prove below that the determinant of this matrix 
is zero when the size of the 0 block is larger than the B 
block (i.e. N > M). 

First, we define a set of N + M vectors e* = 
(0, 0, ..., 1 ,..., 0), such that only the ith component of 


the vector ej is nonzero, while i = 1,2 ,...,N + M. 
Here, these vectors span a N + M-dimensional linear 
space. This linear space is the direct sum of two sub¬ 
space Vj ®V- 2 : where V\ is spanned by the vectors e, with 
1 < i < N and Vj by those with N + l<i<N + M, 
and it is easy to check that V\ and V 2 are orthogonal to 
each other. 

Because the upper-left block of the matrix only con¬ 
tains zeros, etNiej = 0 for i < N and j < N. It 
implies that for any vector ej with j < N, the vector 
Mej is orthogonal to any vectors in the subspace V\. In 
other words, Mej is a vector in the subspace V- 2 - For 
j = 1, 2,..., IV, Mej generates N vectors in the space of 
V 2 . If N > M, some of these N vectors must be linearly 
dependent, because the space V 2 can only have M inde¬ 
pendent vectors. Therefore, we can construct (at least) 
one zero vector using these N vectors A4ej: 

N 

= 0 (G2) 

0 =1 

where a J are N numbers. This implies immediately 
that the M. matrix has (at least) one null vector v = 

2—*j=l a j e j 

N 

Mv = ^ Nlcijej = 0 (G3) 

i=i 

Having a zero eigenvalue implies that the determinant of 
M must be zero. 


Appendix H: isostatic condition, elasticity and local 
vertex-face correspondence 

In this section, we reveal the connection between the 
isostatic condition and the criterion for the existence of 
a local vertex-face correspondence. 

The idea of isostaticity plays an important role in the 
study of mechanics stability. It comes from the counting 
argument developed by Maxwell 61 . If we construct an 
elastic system by connecting beads with rigid rods (and 
allow the rods to rotate freely around each joint), the 
rigidity of the system can be determined by comparing 
the total number of constrains and the total number of 
degrees of freedom. In 2D, the total number of degrees 
here is 2 times the number of beads, because each bead 
has two degrees of freedom in 2D, while the number of 
contains is the number of rods, since each rod enforces 
one constrain by fixing the distance between two beads. 
If we consider such a system as a graph (i.e. beads as 
vertices and rods as edges), the number of degrees of free¬ 
dom is 2 N v , while the number of constrain is N e . The 
isostatic condition requires these two numbers to coin¬ 
cide. If all the constrains are independent (i.e. no re¬ 
dundancy), this condition represents the verge of a me¬ 
chanical stability (i.e. a phase transition point). If we 
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add/remove one edge (rod) to the system, the system 
becomes stable/floppy. The rigorous formula for the iso¬ 
static condition in 2D is 

2N V =N e + 3 (HI) 

Here, a extra number 3 is introduced to the r.h.s. to rep¬ 
resent the trivial global degrees of freedom (two trans¬ 
lations and one rotation), which will alway arise. Using 
the Euler characteristic (N v — N e + Nf = 2 — 2 g), we can 
rewrite the isostatic condition as 

N v = N f + 1 + 2 g, (H2) 

where g is the genus of the underlying manifold. In the 
thermal dynamic limit (N v — > oo and Nf —* oo), we can 
ignore the finite part 1 + 2 g and therefore, the condition 
coincides with the our requirement of N v = Nf. 

In addition, for an isostatic system, to ensure that 
there is no redundant constrains, one shall require that 


for any subsystem (subgraph), the total number of de¬ 
grees of freedom always exceeds (or equal to) the number 
of constrains (plus three) 

2Ad > Ad + 3 (H3) 

where Ad and N' e are number of vertices and edges in a 
subgraph, while 3 on the r.h.s. comes from global trans¬ 
lations and rotations. If a subgraph has the topology of 
a disk (i.e. the Euler characteristic is N v — N e + Nf = 1 ), 
we can rewrite the condition as 

N' v > Nf + 2 (H4) 

which is very similar to but slightly stronger than our 
criterion of local vertex-face correspondence (Ad > Nf). 

Because our criterion is slightly weaker, some of the 
lattices/graph that are not isostatic can still be used to 
construct a discretized Chern-Simons gauge theory, e.g., 
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